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FOREWORD 


Lockheed Missiles & Space Company, Inc., Huntsville Research 
& Engineering Center, submits this final report to NASA-Marshall 
Space Flight Center as fulfillment of Contract NAS8-35987. The 
NASA technical monitor for this contract was Klaus W. Gross, EL24, 
of the Systems Analysis and Integration Laboratory. 
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1. INTRODUCTION 


During the past few years the requirement for more accurate and econom- 
ical computational methods for predicting fluid flows in practical engineer- 
ing configurations has become increasingly important. Full-scale testing 
or model testing is often expensive and in many cases inconclusive, necessi- 
tating heavier reliance on computational procedures. It is also important 
to have adequate methodologies and computer codes available to help in the 
interpretation of experimental data and to aid in fluid flow design prob- 
lems . Many computer codes and associated methodologies have been generated 
to provide solutions to complex geometric flowfield configurations. The 
treatment of fluid flows by these finite differencing methods has advanced 
tremendously in recent years due to production of advanced machines such as 
the CDC CYBER 205 and the CRAY. 

The current computational method for determining rocket thrust chamber 
performance is the JANNAF procedure. This utilizes the TDK and BLIMP com- 
puter codes. Through the advances in computer technology, new codes have 
been developed which solve the Navier-Stokes equation with finite difference 
or finite element solution techniques. Since this technology represents 
flow fields in general, solutions may be obtained for many conditions that 
exceed the existing TDK/ BLIMP capability. 

This report presents the results of a contract effort by Lockheed- 
Huntsville to compute the viscous, axisymmetric flow in the SSME thrust 
chamber. The Lockheed-Huntsville developed General Interpolants Method 
(GIM) code was used for the calculation on the CYBER 205 computer at the 
NASA- Langley center. The results of the study show that Navier-Stokes codes 
can indeed by used for these flows to study trends, viscous effects, and to 
determine flow patterns. Limitations on our method for absolute accuracy of 
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perf oraance numbers have been noted and are explained in this report . The 
conclusions are summarized as follows: The Navier-Stokes codes that we are 

familiar with including the GIM code can be used for qualitative study of 
viscous nozzle flows, but further research and development is needed before 
they can be used as production tools for nozzle performance calculations. 

Section 2 describes the GIM formulation, numerical scheme and computer 
code that was used in the study. Section 3 discusses the actual SSME nozzle 
computation, showing grid plots, flow contours and flow parameter plots. 

The computer system and run times/costs are detailed in Section 4. Conclu- 
sions and recommendations for further analysis are given in Section 5. 
Appendix A presents the detailed flowfield solution in computer printout 
form including a description of the mesh and identification of flow 
variables . 
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2. GIM CODE DESCRIPTION 

2.1 THE GIM METHODOLOGY 


The General Interpolants Method (GIM) code (Refs. 1 through 7) was 
developed to analyze complex flow fields which defy solution by simple 
methods. The code uses numerical difference techniques to solve the full 
three-dimensional time-averaged Navier-Stokes equations in arbitrary geo- 
metric domains. The equations are cast in strong conservation law form and 
written in an orthogonal Cartesian coordinate system. Included are a con- 
tinuity equation for global mass conservation, three components of momentum 
conservation, total energy conservation and equations for conservation of 
individual species. Pressure is related to the conservation variables 
through a gas law for arbitrary mixtures. A generalized geometry package is 
used to model the flow domain, generate the numerical grid of discrete 
points and to compute the local transformation metrics. Computation is done 
in physical space by explicit finite difference operators. The GIM approach 
essentially combines the finite element geometric point of departure with 
finite difference explicit computation analogs. This provides a capability 
which takes advantage of the geometric flexibility of an element description 
and the computation speed of difference representations. 


The numerical analogs of the differential equations are derived by 
representing each flow variable with general interpolant functions. Con- 
sidering a general quadrilateral element with four nodal points, the differ- 
ential equations are integrated over this element at each nodal point using 
Gauss quadrature. See Fig. A-l in Appendix A for an example of a typical 
two-dimensional grid element. The point of departure then requires that a 
weighted integral of interpolants be zero over the flow domain. By choosing 
the weight functions to be the interpolants themselves, the GIM formulation 
produces identically the classical implicit finite element 
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discrete equations. These forms are not used in the GIM code due to their 
fully implicit nature and inherent inefficiencies. Rather, the weight func- 
tions are chosen to be orthogonal to the interpolant functions which pro- 
duces explicit finite difference type discrete analogs. By appropriate 
choice of constants in the weight functions, the GIM becomes analogous to 
standard finite difference schemes such as centered, backward, forward, 
windward and multi-step predictor-corrector schemes. The GIM analogs, how- 
ever, are automatically produced for arbitrary geometry flow domains and 
hence are a general point of departure and provide flexibility in the choice 
of differencing schemes. The following discussion summarizes the GIM formu- 
lation. 

2.2 GEOMETRIC TREATMENT 

The GIM formulation is not a finite element method in the classical 
sense. Instead, finite difference methods are used exclusively but with the 
difference equations written in general non-orthogonal curvilinear coordi- 
nates. Transformations are used to transform the physical planes into 
regions of unit cubes. The mesh is generated on this unit cube and the 
local metric coefficients generated. Each region of the flow domain is 
likewise transformed and then blended via the finite element formulation to 
form the full flow domain. In order to treat "completely— arbitrary" geo- 
metric domains, different transformation functions can be employed in dif- 
ferent regions. The blended domain is then transformed back to physical 
space and the Cartesian set of equations is solved for the full region. The 
geometry part of the problem is thus treated much like a finite element 
technique, while integration of the equations is done with finite difference 
analogs. 

An important concept of the theory of finite elements is the "discon- 
nectness property. This means that a domain can be divided into a finite 
number of pieces called elements and the approximation to the functional 
distribution over each element can be studied independently. Thus the 
approximating functions for each element completely define the behavior of 
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the function profile within the element without consideration of its ulti- 
mate location in the full model. After each element is defined, the com- 
plete discrete model of the body is obtained by “assembling the system." 
The assembly is performed by means of mathematical blending of each sub- 
domain while maintaining continuity at the junctions. 


The GIM code uses these concepts, borrowed from finite element theory, 
to obtain discrete finite difference models of the Navier-Stokes equations 
in arbitrary geometric domains. The development is done in local curvi- 
linear intrinsic coordinates based on the following concepts: 


• Analytical regions such as rectangles, spheres, cylinders, 
hexahedrals, etc., have intrinsic or natural coordinates. 

e Complex regions can be subdivided into a number of smaller re- 
gions which can be described by analytic functions. The de- 
generate case is to subdivide small enough to use very small 
straight-line segments. 

• Curvilinear intrinsic coordinate systems result in constant 
coordinate lines throughout a simply connected, bounded domain 
in Euclidean space. 

• The intersection of the lines of constant coordinates produce 
nodal points evenly spaced in the domain. 

e Curvilinear intrinsic coordinate systems can be produced by a 
univalent mapping of a unit cube onto the simply connected 
bounded domain. 


Thus, if a transformation can be found which will map a unit cube uni- 
valently onto a general analytical domain, then any complex region can be 
piecewise transformed and blended using general interpolants . 


2.3 GIM ELLIPTIC MODULE 

The partial differential equations solved by the GIM elliptic module 
are the Navier-Stokes written in three-dimensional Conservation law form for 
a Cartesian coordinate system. These are shown in Fig. 1. 
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Fig. 1 Navier-Stokes Equations 
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In the equations, T represent the viscous stress terms, q the heat con- 
duction and R the mass transfer contributions. The viscosity coefficients, 
M, the thermal conductivity, k, and the binary diffusion coefficient ,£) 
are obtained by various closure laws for viscous flows. For inviscid, 
hyperbolic flows with shock waves, these stress terms are also used, but 
with artificial viscosity coefficients, to facilitate capture of strong 
shocks. The GIM formulation has built-in formulas for computing the arti- 
ficial component of the coefficients. These are termed Numerical Diffusion 
Cancellation (NDC) coefficients. The purpose of these NDC coefficients is 
to cancel low order truncation error terms which arise in the numerics. 

The GIM code solves these equations using explicit finite difference 
methods. This is done with a set of general interpolation functions iOl-^, 
such that 


f( WV = 2 VvW f j 

where f^ are the flow variables at the corner points of the element. The 
simplest form for the are the trivariable Lagrange interpolants . The 
theory itself does not restrict the I ^ to be linear, but in the present 
code we currently have only the trilinear interpolants. Any flow gradient 
can then be simply computed as 


3f = 
9n i 



91 . 

— - f . 
*i\ ± j 


The code uses these concepts to produce a discrete analog of the 
Navier-Stokes equations for a single analytical region. The point of 
departure is the requirement that 


/ 


W $ dV 


0 
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where W is an arbitrary weight function and 0 is the differential equation 


$ 


9t 9x 9y 9z 


The general interpolants I are used to approximate each of the flow vari- 
ables U,E,F,G and are used to obtain a discrete analog; 

[A. .] U. + [B, J E. + [C, .] F. + [D. .] G. = 0 

ij J +J J ij J ij J 

where Uj,E^,Fj,G. are the flow variables at the node points of an element. 
The coefficient matrices are geometrically dependent only. It is this 
general finite difference expression that is solved in the GIM code. 

2.4 GIM PARABOLIZED MODULE 


The full three-dimensional Navier-Stokes equations just discussed are 
elliptic in character and require either time-dependent or relaxation/ 
iteration schemes to integrate the complete spatial flow field simultane- 
ously. This can require large amounts of computer storage and relatively 
long run times. If the physical problem is of the nature that such a solu- 
tion is required, then one has no option other than to "bite the bullet" and 
use it. 

For situations in which a region of the flow is inviscid and entirely 
supersonic, a spatial hyperbolic marching algorithm would be efficient. 

There are also many viscous problems of interest in which a parabolic march- 
ing solution would be acceptable. The most versatile and efficient use of a 
GIM code is to combine all three options into a single user-oriented code. 

In this manner the elliptic algorithms can be called upon to compute in 
elliptic flow domains up to a hyperbolic downstream boundary. Automatic 
and/or user supplied conditions can then be used to switch to the forward 
marching routines to compute the remaining flow domain. 
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The basic idea of the GIM parabolic technique (Ref. 6) is to combine 
the classical parabolic marching approach with a "quasi-time" relaxation. 
This will be termed the GIM parabolized module. The parabolic -march pro- 
cedure greatly reduces the amount of computer storage compared to a fully 
elliptic field. The time relaxation form of the equations eliminates the 
decode ambiguity associated with the parabolic pressure problem and allows 
velocity boundary conditions at solid walls to be treated. The equations 
used in the parabolized module are the time-averaged full Navier-Stokes , but 
with all second order terms dropped in a quasi-marching coordinate. Another 
way to view the equations is to take the parabolized Navier-Stokes and add 
back pseudo time derivatives. The solution procedure, as any parabolic 
marcher, thus allows no streamwise diffusion effects. The solution is 
assumed known at upstream data planes, 1, 2, ... K-l, and the solution is 
sought at plane K with no knowledge of plane K+l. Quasi-time relaxation is 
used to obtain the solution at plane K in terms of the (converged) solution 
at a number of upstream data planes. Second order backward streamwise dif- 
ferences are used to prohibit downstream feedback. Thus the parabolized 
algorithm is not a classical space marching scheme, and is also not a time- 
dependent elliptic method. It is somewhat of a hybrid technique which com- 
bines the better features of two approaches and eliminates some of the bad 
ones . 


The procedure for marching a solution is outlined as follows: 


• A starting plane of data and gas properties are provided at 
k=T . 

• The first calculated plane k=2 is obtained from a first order 
finite difference equation. 

• An initial guess is provided at t=0 for all flow variables at 
plane k=3 . 

• The Predictor-Corrector is iterated to steady state 

U + 0 

which provides the parabolized solution at plane k. 
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• The conservation variables are decoded and the plane k is set 
k+1. The marching procedure is then repeated. 

2.5 FINITE DIFFERENCE SCHEME 

The finite difference scheme used in GIM can best be described as three 
point "backward in x" and "MacCormack in the y-z plane." The system is en- 
tirely explicit in that all nodes are initialized at pseudo time t=0 and 
relaxed to steady state using known or previous data. The grid and differ- 
ence equations used for illustration are Cartesian-like for simplicity; the 
code, of course, has the geometric generality described earlier. All data 
are known for planes k-2, k-1 and the solution is to be marched to plane k. 
The following two-step predictor-corrector sequence is used: 

Predictor 


U n+1 

i.j.k 


I? . . 


At 

Ax 


JX 


E — j? 


1 At 

2 Ax 


* * 

E. . . , - E 
i,j,k-l i,j,k-2 


J 


At 

Ay 


r F n _ F n ~i _ At 

c n - r n 

i+1, j,k i,j,kj Az 

^i,j,k 


Corrector 


U 


n+1 

i, 


1 TJ n ^.n+1 3_ Ait ^n+1 

2 i,j,k i,k,k 2 Ax i,j,k 




1 &t r * 

2 Ax i,j,k-l 


* 

E i,k,k-2 


At 

Ay 



F 


n+1 1 

i-1, j,k 


At c;n+l _ p,n+l 
Az i,j,k i,j-l,k 


The "hat" (~) indicates provisional values from the predictor step and the 
superscript (*) denotes converged values at upstream planes k-1, k-2. 
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Stability of this explicit scheme is obtained by obeying the classical CFL 
constraint of explicit methods. An artificial viscosity is added to the 
equations for stability when capturing a strong shock. 

2.6 SUMMARY OF THE GENERAL INTERPOLANTS METHOD 






m 

© 


Numerical Solution of Conservation Laws 

Based on Time-Average Navier-Stokes Equations 

Elliptic or Parabolized Solution Algorithms 

Geometrically Arbitrary with Finite Element-Type Treatment 

Grid Generation Using Trivariate Blending Function Interpola- 
tion 

Geometry Input by Interactive Graphics 

Finite Difference Nodal Analog Solution Algorithms 

Specific Scheme Selected by User Input 

Boundary Conditions via Finite Difference for Internal or 
External Flows 

Integration by Predictor-Corrector Scheme 
Shock Capturing Approach for Supersonic Flow 
Ideal and Real Gases 
Inviscid or Viscous Flows 

Laminar Viscosity or Algebraic Eddy Viscosity Turbulence Models 
Two or Three-Dimensional Capability 

Code Operational on Univac 1108, CDC 7600, and Vectorized for 
CDC CYBER 203 and 205. 


2.7 STRUCTURE OF THE CODE 

The GIM code is divided into four modules: (1) mesh generation; (2) 

nodal analog assembly; (3) unsteady integration; and (4) data display. The 
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mesh generation module accepts boundary geometry data, curve or line formula 
flags, and number of cuts in each coordinate direction. A set of general 
curvilinear coordinate maps is then used to subdivide each region into 
finite elements. Each region which is input is processed and then blended. 
The output is a set of coordinates for each element along with the element 
coefficient matrices. The nodal analog assembly module takes the mesh data 
from a stored external file and performs, via quasi-variational procedure, 
the assembly of the element equations into the full domain equations. At 
this point, the dynamic storage allocation is set up so that the unsteady 
integration module can integrate with virtually unlimited problem size. 

The unsteady integration module performs the actual computation of the 
flow by employing the boundary conditions selected by the user. The nodal 
analog at this point is arbitrary and any one of a number of schemes can be 
selected depending on the problem being analyzed. The solution is marched, 
forward in time for a specified number of steps or until a steady state is 
reached. The data display module reads the solution profiles from external 
storage (drum/tape) and prints, plots and maps the flow parameters. 

The GIM code is currently operational on three computer systems: 

Univac 1108/1110 (2-D Scalar, Small Inviscid Problems) 

CDC 7600 (2-D Scalar, Small Viscous Problems) 

CDC CYBER 203 and 205 (3-D Vectorized, Viscous Problems) 

Representative run times for the elliptic code are given in Chart 1. 


System 

2-D 

3-D 

Example 

10,000 Node, 500 Steps, 3-D 

Univac 1108 

3.7 x 10" 3 

5.9 x 10~ 3 

8.1 hours 

CDC 7600 

3.6 x 10 ^ 

5.8 x 10 ^ 

48 minutes 

CYBER 203 

2.2 x 10~ 5 

3.3 x 10~ 5 

3 minutes 

CYBER 205 

5.5 x 10’ 6 

8.3 x 10“ 6 

1 minute 


Chart 1 CP Run Time for GIM Code (Seconds for Node per Iteration) 
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3. SSME THRUST CHAMBER PERFORMANCE CALCULATION 
3.1 GOVERNING EQUATIONS 


The symmetry of the Space Shuttle Main. Engine (SSME) nozzle flow about 
its axis can be used for better computational economy and greater resolution 
by using an axisymmetric numerical approximation. The complicated three- 
dimensional Navier-Stokes equations in Section 2 can then be simplified to 
the following axisymmetric form (Ref. 8): 


where 


(Y 

■~r + + - -- (rE^) + H W = 0 

3t 3x r 3r 2 


U = (Ojpu^ pu 2 ,p<?) 


_ TT a ^ .0(2 .0(3 . ^ 

E7 = u.U - t 6 - t. 0 6 ~ (u.t. . + q.) 

i l ll i2 2 ij i 


H a = [9iT/9r + 2uu 0 /r*] 6 


2, .0(3 


= ^(Bu^/ax. + auVax^ - 

TT = P - X(3u i /ax i + u 2 /r) 
q i = k 3 T/ax i -^2^6.2 

and summation over repeated indices is implied. 

In the equations, p, u, t, T, q, y , and A have the same meaning as 
those in Section 2. (E^, thus, is equivalent to F.) The symbol 6^ or 
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Is the Kronecker delta, the superscript a is the equation index = 1,...,4 
and the subscripts i,j are the spatial coordinate index corresponding to 
x,r, the axial and radial coordinates. 

A variation of the Petr ov-Galerkin explicit finite element scheme for 
axisyrametric viscous compressible flows has been developed and incorporated 
into the General Interpolants Method (GIM) code. This variation preserves 
the second order accuracy in time and space of the regular Petrov-Galerkin 
method even in the radial direction without the restriction of r > > A r. 

More detail is available in Ref. 8. 

3.2 SSME NOZZLE GEOMETRY 

•k 

The SSME thrust chamber has a throat radius, r , of 5.1527 in. and a 
total length of 11.167 ft, of which 10.10 ft is the length of the divergent 
section. Based on the wall contour data points obtained from Klaus W. Gross 
of NASA- Marshall Space Flight Center the nozzle wall from the tangent 
point, x w , downstream can be curve-fitted as: 


r = -32.44705 + [1123.9917 + 102.1846 (x - x )] 1/2 - 0.77040726 (x - x ) 

w 1 w 

where 

^ ^ 

r = r/r , x = x/r , x = 0.23591149 

w 

The wall contour of the nozzle and the coordinate system are shown in Fig. 

2. The throat area detail is shown in Fig. 3. 

To model the flow domain of the SSME nozzle as defined by Figs. 2 and 3 
a numerical grid of discrete points was generated. The choice of the grid 
size is primarily dictated by the solution accuracy, convergence, and 
stability criteria. The final grid used in the computation is shown in Fig. 
4. It can be seen that in the regions close to the wall and around the 
throat the grid spacing is much finer than that in other regions. Figure 5 
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- - - 1/2 
r = -32.44705 + [1123.9917 + 102.1846 (x - x )] 

w 

- 0.77040726 (x - x ) 

w 



r = 8.8034 x = 23.5232 

e e 

x = 0.23591149 

w 


Fig. 2 Axisymmetric SSME Nozzle 
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Fig. 3 SSME Nozzle 
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Fig. 4 Axisymmetric SSME Nozzle Grid 



Fig. 5 Enlarged View of Grid near Throat Region 
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shows the throat region enlarged. The grid spacing change was carried out, 
where practicable, in a manner to minimize any sudden spacing variation. 
These considerations were necessary to obtain better resolution in the 
numerical results. There are 6480 (36 x 180) grid points in the final grid 
layout used in the computation. The numbering scheme for the computational 
mesh is given in Appendix A. 

3.3 INITIAL CONDITIONS 


The steady state solution of the SSME nozzle flow problem was found as 
the limiting approach of the time dependent solution to the governing equa- 
tions noted in Section 3.1. In order to obtain a rapid convergence po the 
steady state solution the initial values of density, velocity, pressure and 
other variables were determined from a one-dimensional isentropic expansion. 

I 

They were derived from the Chemical Equilibrium Composition (CEC) code of 
NASA-Lewis Research Center (Ref. 9) using 100 percent SSME power level 
operating conditions: 


Chamber Gas Stagnation Pressure, psia: 2935.7 

Chamber Gas Stagnation Temperature, R: 6550.2 

Oxygen/Hydrogen Ratio: 6.0549 

Nozzle Expansion Ratio: 77.5 

Liquid Enthalpy, cal/mol: -1837.660 

Liquid Enthalpy, cal/mol: -2884.385 


The results from the CEC calculation, assuming equilibrium composition dur- 
ing expansion through the nozzle, are shown graphically in Figs. 6 through 
9. They are also shown in tabular form in Table 1. This table was used in 
the initialization for the GIM code computation using a table look— up and 
linear interpolation procedure to assign an initial value to each grid point 
according to its axial location. It is emphasized that the isentropic ex- 
pansion exponent (Y) and the molecular weight (MW) of the oxygen/hydpogen 

* 

mixture were varied as a function of x/ r as shown in Fig. 9 during all 
the subsequent numerical calculations. 
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Fig. 7 Total Velocity Initial Values from CEC Solution 


20 


LOCKHEED-HUNTSVILLE RESEARCH & ENGINEERING CENTER 






Weight 


LMSC-HREC TR D951729 



O Molecular Weight 


Fig. 9 Molecular Weight and Gamma Initial Values from CEC Solution 
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Table 1 FLOWFIELD INITIALIZATION DATA FROM CEC 
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3.4 BOUNDARY CONDITIONS 

The SSME flow field was computed using the elliptic module of the GIM 
code. The boundary conditions consist of a subsonic inflow plane with total 
pressure and total temperature prescribed, a symmetry axis condition at the 
nozzle centerline, and a supersonic (one-sided) outflow condition at the 
nozzle exit to simulate a near-vacuum expansion. At the wall, in addition 
to the viscous no-slip condition, the wall temperature boundary condition 
was determined from the CEC analysis. The input equilibrium chemistry wall 
temperature profile from the CEC analysis is shown in Fig. 10. Figure 11 
shows schematically the boundary conditions used in the computations. 

3.5 TURBULENT VISCOSITY MODEL 

Based on the initial values of the density, velocity, and viscosity de- 
rived from the CEC calculation the Reynolds number of the flow, Re^, based 

7 u 

on the SSME nozzle throat diameter, is in the order of 10 . The flow is 
thus expected to be turbulent. The turbulent model used in the computation 
was based on Baldwin-Lomax model (Ref. 10) modified by P.D. Thomas for the 
nozzle wall boundary layer. The modification mainly concerned the method 
used to compute the vorticity function (Ref. 11). The viscosity used in the 
computation is the sura of the laminar viscosity from the GIM code input and 
the turbulent viscosity calculated from the Baldwin-Lomax model. In this 
model, the eddy viscosity length scale is determined from the spatial dis- 
tribution of vorticity. Therefore, the turbulent viscosity in the region 
far away from the nozzle wall boundary layer is nearly null, and the total 
viscosity in that region will be essentially laminar. This pertinent fact 
will be demonstrated by the following brief description of the model. 

The turbulence model used here is a two-layer algebraic eddy viscosity 
model in which the inner viscosity is defined by the Prandtl-Van Driest 
formula : 
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Subsonic Inflow Condition 
with Prescribed Total Pres- 
sure and Total Temperature 


Fig. 11 Boundary Conditions 
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(u ) = p * 2 M, 

inner 


where, ^{denotes the normal distance from the wall, 

A - 0.4 ^[1 - exp(-^ + /26)], 

=>efi Re p w | to | w /p w )* 

and,^J^ is the smallest value of^j^at which the inner and outer viscos- 
ities are equal. The subscript w denotes conditions at the wall, and |uj| i 
the magnitude of the vorticity. The corresponding outer viscosity is given 
as: 


(y t ) = 0.0168 C . P F • F k fe£>, J. >J 

outer v e 


where, C C p is a constant, and 


„ ) wk. nr-^m m 

F = ' * 


1^4 f m* 


u < 2 f 

m m 

u > 2 f 

m ~ m 


\<d> = 


2 

i 

m 



The quantities-^ and f are determined from the function 
mm 

fJcf) = >cf |w| {1 “ exp(-^/ + /26) ] . 


The following values of the constants appearing in the foregoing relations 
were used in the computation: 
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C = 2.08 

cp 

C k = 1 ‘° 

C wk = °* 25, 

The viscosity profile at the exit plane as a function of the distance from 
the nozzle wall is shown in Fig. 12. At other cross-sectional planes the 
viscosity profile is similar. The laminar viscosity (h^ am ) used in the 
computation was 0.74 x 10 ^ lb/ft-sec. It can be seen from Fig. 12 that 
in the region from the axis to the half way point from the wall the vor- 
ticity is negligible and the viscosity is laminar, and only in the boundary 
layer and in the region from the boundary layer edge to one half the nozzle 
wall radius does the turbulent viscosity become dominant. 

3.6 CALCULATED FLOW FIELDS 

Using previously described boundary conditions and the turbulent vis- 
cosity model, the GIM code was run from the initial conditions described in 
Section 3.3 until a steady state solution was reached. The solution for the 
case with equilibrium chemistry wall temperature profile assumption is shown 
on the figures appearing in the following pages. They are: 

Figure No. Description 


13 

Mach 

Number Contours 

14 

Pressure Contours 

15 

Mass 

Flow through the Nozzle 

16 

Axial 

. Mach Number Variation 

17 

Wall 

Pressure Variation 

18 

Wall 

Density Variation 

19 

Exit 

Plane Pressure Profile 


It can be seen from the Mach number and pressure contours (Figs. 13 and 
14) that the solution in the subsonic/transonic region is very good. In the 
supersonic region the solution is good except in the neighborhood of the 
axis. The contour plots clearly show the development of a thin boundary 
layer adjacent to the nozzle wall. The boundary layer thickness is about 
one grid width. The oscillation of the Mach number and pressure near the 

28 


LOCKHEED-HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 


1.0 


3 

u 



10 5 10" 4 10" 3 10" 2 10" 1 10° 


^Total (lbm/ft-sec) 
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Fig. 12 Viscosity Profile at Nozzle Exit Plane 
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axis can also be seen in Figs. 16 and 19. We believe that the problem near 
the axis is due to the numerics and the coarse grid used there. The mass flow 
through the nozzle is nearly constant with a small variation in the subsonic 
region (Fig. 15). The mass flow rate at the throat is 978 lb/sec. The varia- 
tion in mass conservation from the chamber to the exit plane is approximately 
+0.8 percent. Figure 17 shows the wall pressure variation. The pressure at 

the nozzle wall lip is approximately 1020 lb/ft (7.1 psia) which is some- 

2 

what higher than 780 lb/ft (5.42 psia) as calculated using the Method of 
Characteristics solution. The thrust, the sum of the flow rate "dotted" into 
the velocity and the static pressure times the throat incremental area each 
term evaluated at the throat plane plus the axial component of the wall static 
pressure times the wall incremental area evaluated from the throat to the exit 
plane, was calculated as 512000 lbf and the specific impulse is 523 sec. Due 
to the high wall pressure the thrust and specific impulse numbers just re- 
ported are unreasonably high. These are not to be construed as the final SSME 
performance numbers. 

The problem with the 2-D axisymmetric analysis that has been performed is 
the calculation of the axis conditions. The axis calculation, using 
predictor-corrector techniques, is inaccurate at small values of r. This 
leads to truncation errors which require numerical damping to control, result- 
ing in lower than expected Mach number at the axis and larger than expected 
wall pressure at the nozzle lip. The axisymmetric treatment of the axis in 
the Navier-Stokes codes that we are aware of in the open literature clearly 
require further work before accurate thrust chamber performance numbers can be 
generated. 

A detailed computer printout for the steady state flow field and a de- 
scription of the printout are given in Appendix A. 

The RFP for this effort requested that the heat transfer to the SSME 
nozzle wall be calculated. The GIM code does not presently calculate the heat 
transfer coefficient. However, all the elements required to calculate a con- 
vective heat transfer coefficient are given in the computer printout found in 
Appendix A. 
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Fig. 15 Mass Flow Through the Nozzle with Equilibrium 
Chemistry Wall Temperature Profile Assumption 
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Fig. 16 Mach Number Along the Nozzle Axis with Equilibrium Chemistry 
Wall Temperature Profile Assumption 
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Wall Pressure (lb/ft' 
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Nozzle Length (ft) 


Fig. 18 Wall Density Variation with Equilibrium Chemistry 
Wall Temperature Profile Assumption 
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Fig. 19 Exit Plane Pressure Profile with Equilibrium 
Chemistry Wall Temperature Profile Assumption 
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4. COMPUTER REQUIREMENTS 

4.1 COMPUTER POWER REQUIRED FOR VISCOUS NOZZLE CALCULATIONS 

There are two basic concepts in modern digital computer architecture: 
scalar and vector. A scalar machine performs operations (adds, multiplies) 
sequentially, one after the other. A vector machine performs ’’parts" of 
multiple operations all at the same time. On the CYBER 205 for example, 
65,000 multiples can be going on at the same time. This is obviously many 
times faster than a scalar operation. All Class 6 supercomputers use the 
vector concept while the older Univac 1100 and CDC 7600 machines are based 
on scalar architecture. 

The GIM code elliptic version used in this analysis requires the use of 
a large vector processor. The NASA-Langley CYBER 205 with its 16 million 
word memory is one of the most powerful computers available. This is the 
type of machine that will be required for the solution of viscous nozzle 
problems using Navier-Stokes solution techniques. The requirement of fine 
mesh near the nozzle wall to resolve the boundary layer and fine mesh in the 
transonic region led to the selection of a large vector machine such as the 
CYBER 205 as one of the few possibilities for the solution of viscous in- 
ternal flow problems. Scalar codes can be used to solve this type of prob- 
lem but are not cost effective for viscous nozzle calculation. 

4.2 COMPUTER COSTS FOR THE GIM CODE SSME VISCOUS NOZZLE CALCULATION 


The solution reported in Section 3 of this report required 20,000 iter- 
ations to reach steady state. The performance of the CYBER 205 is measured 
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in Computer Resource Units (CRUs) which is approximately one second of exe- 
cution time. The SSME nozzle calculation required 0.2 CRUs per iteration 
tor a total of 4000 CRUs per case. The rate at the NASA-Langley Computer 
Center is approximately $0.75 per CRU . Allowing a few CRUs contingency and 
including the computer time required to set up the computational mesh, each 
case is estimated to cost approximately $3500 to compute. 
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5. CONCLUSIONS AND SUMMARY 

5 .1 THE STATE OF THE ART OF VISCOUS NOZZLE CALCULATIONS 

Performance calculations in large area ratio nozzles remains as 
one of the challenges of CFD Navier-Stokes codes. A number of critical 
items remain unsolved before the research codes of today can be used to 
design nozzles of the future. Among these are the following: 

1. Algorithms for treating the ’’axis" singularity in the axisym- 
metric equations. 

2. Turbulence models for treating the wall boundary layer in tran- 
sonic and supersonic flow. 

3. Subsonic boundary condition treatment at the "inflow" boundary 
to allow numerical waves to pass through it. 

4. Numerically stable algorithms which do not require large arti- 
ficial damping. 

5. Algorithms for switching between elliptic solution in transonic 
flow to parabolic solution in supersonic flow. 

6. Efficient models for treating chemistry (equilibrium and finite- 
rate) in viscous nozzle flow. 

The axisymmetric equations for viscous flow contain 1/r terms which 
become indefinite at the axis (r = 0). Most algorithms in use today, 
MacCormack, Euler, etc., introduce false losses into the solution in an 
attempt to resolve the axis flow. The problem gets worse as the grid is 
refined (Ar 0) since l/(r Ar) appears in the finite difference equations. 
One solution is to treat the full three-dimensional equations which avoids 
the singularity. This is however, very inefficient in axisymmetric nozzles 
and requires many times more computer resources than the 2-D, axisymmetric 
equations. Basic research into numerical algorithms for treating the axis 
singularity is needed before actual performance calculations can be made. 
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Computation of viscous flow at high Reynolds number requires a specifi- 
cation of the turbulent viscosity coefficient. Many ad hoc models currently 
exist for this specification. Many arguments also exist as to which model 
is "best" for nozzle flows. The two-equation k-e turbulent kinetic energy 
models seem to be the most exact. However, they require large amounts of 
computer time and the wall boundary conditions for these equations are 
unknown. Algebraic models, termed "eddy viscosity" are the most practical. 
The ones which use vorticity as the governing parameter have worked the best 
to date. Their accuracy in separated boundary layers remains questionable. 
Perhaps the only real way to ever treat turbulence is to solve the full 
Navier-Stokes equations including all the "fluctuating" terms rather than 
the Reynolds-averaged equations. Very large supercomputers, with very long 
run times will be required. 

The inflow boundary condition in an internal, subsonic flow is a crit- 
ical item. Mathematically, four boundary conditions are required in axisym- 
metric flow; three at the inflow and one at the outflow. Specifying all 
four at the inflow boundary results in an ill-posed problem with numerical 
waves stacking-up against the boundary. In the viscous, axisymmetric nozzle 
flow, total pressure, total temperature and flow direction are good choices 
at the inflow, with static pressure, or density (mass flow) at the transonic 
(throat) boundary. At the supersonic nozzle exit, no boundary condition is 
required (except perhaps in the subsonic boundary layer region near the 
lip). The current state-of-the-art in CFD has not checked out this proper 
treatment. 

Artificial damping is currently required in all CFD algorithms. This 
"numerical viscosity" is present in all schemes whether it is added by the 
user directly or just simply occurs due to truncation of the infinite series 
representation. If enough grid points could be used in a high Reynolds num- 
ber boundary layer, then large damping would not be needed. However, this 
is very impractical in terms of computer time. An algorithm which treats 
the "thin" boundary layer without large numbers of grid points is needed. 
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Proper and efficient treatment of the full nozzle is best done with an 
elliptic solution in the subsonic/ transonic region and a parabolized Navier- 
Stokes solution in the supersonic region. Currently available codes require 
the user to switch the algorithms manually. Certain smoothing logic is also 
needed to blend the solutions. An automated procedure and code could be 
constructed which would perform this switching in a production code. 


An accurate chemistry treatment is also essential if performance num- 
bers are to be predicted. Currently available equilibrium models range from 
table look-up of previously obtained data as we have done in this solution 
to in-line solution of the reaction equations. A combination of equilibrium 
and frozen chemistry is usually required. Solution of nozzle flows with 
non-equilibrium chemistry can be made with existing codes, but only for 
benchmark cases. Nozzle design using these chemistry codes is still not 
practical. 

Current CFD Navier— Stokes codes are producing some amazing results for 
a large variety of problems. The future of such codes for nozzle design is 
good, but much research remains to be done. The areas discussed above 
represent a few of these needed research areas. 

5.2 COMPUTER SYSTEM REQUIREMENTS 

The advent of the Class 6 Supercomputer has revolutionized the field of 
computational mechanics. The CFD codes in use today are barely practical on 
the CRAY and CYBER 205 machines. The NASA-MSFC Univac computer cannot be 
used for running these codes to solve two— or three-dimensional viscous 
problems. One of the major problems that we have encountered in SSME flow 
calculations is the computer time purchase situation. Because MSFC does not 
have a supercomputer, our codes must be run at Langley, Ames, or even com- 
mercially to Control Data Corporation. This is very expensive since MSFC 
must send "money" of some kind to the other NASA centers. We strongly 
recommend that MSFC get a supercomputer as soon as possible. We also 
strongly recommend the CDC CYBER 205 and not the CRAY 1. 
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The very large memory of the 205 and the extra CPU speed make it by far the 
best on the market today. 

This most important item of computer time must be solved before MSFC or 
its contractors can economically use CFD codes to calculate the performance 
of advanced propulsion systems. 
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Appendix A 

PRINTOUT AND DESCRIPTION 
OF STEADY STATE FLOW FIELD 
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Appendix A 


A printout of the steady state flow field at iteration 20000 is 
listed on the following pages. The nomenclature for the printout is: 

NODE nodal point number (see Fig. A-l) 

X,Y Cartesian coordinates of the NODE (Y = r), ft 

3 

RHO static density, lbm/ft 

u,v velocity components in X,Y directions, ft/sec 

SOS local sonic velocity, ft/sec 

2 2 

E total energy, ft /sec 

2 

P static pressure, lbf/ft 

IB node point type (see Table A-l) 

GAM local value of the specific heat ratio 

CV constant-volume specific heat, Btu/lbra-R 

THXY flow angle in X-Y plane, deg 

THXZ flow angle in X-Z plane, deg 

MU local viscosity, lbm/ft-sec 

T temperature, R 

M Mach number 
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.116541116+06 

• 1 

.235994745+06 
. 1 01290465+01 
•2375 4277 5+ C6 
. 1 0C7542 ..5 + 01 
.239 099455 + 06 
. 102361895+01 
.239989645+06 
.1 C0345595+01 
*239085085+06 
.100410655*01 
.238558575+06 
.1 00657275+01 
.237608575+06 
.1 01037675+ul 
•236411175+06 
.101526225+01 
. 235 0893 15 + 06 
. 1 72 594376*0 1 
.233352095+06 

• 1- 27°. 6526+01 
. 231621445+06 
.17354 35 5E+J1 
.22°54'’ Of. 6 + 0 6 

• 1044-57 05 + 0,1 
.2 27291165+06 
.1 0-5331725 + 01 
.22480825+06 
. 1 063255 55+ .1 
•222242155+06 
. 1 074 *4025+01 



2 r 22 3 156 5+ ’3 



A-15 


ITERATION 2'fsj 

2 T I M r 

= . 2*767? C ? E 

-06 T I 

‘‘5 = . 3 4 6 7 3 5 9 

r «, 1 2 




NODE 

Y 

Y 

RHO 

! J 

V 

SOS 

5 

D 

13 



C-AM 

CV 

T J X Y 

thxz 

u u 

T 

M 


12 04 

.3245228 7E-.1 

.26 '.36 4 32 5 + 

.316626=75+ u 

. 84 = 27 '38 r + 74 

.22 ■'97939E+ 03 

.5 06451 795+ • 4 

. l6375S94E + -'9 

•21=486085+06 

3 



• 1 1 5 G v 3 7 7 5 + ■. 1 

.=52351775+ . 

• 23 228. 17 = 5 + 71 

# 1 

,74»' iOf- 14 

.62 3 3564 15* '4 

. 1 0=542285+01 


12 05 

.8 J54613 05-01 

• 2 7 1 4 1 7 .0 4 E ♦ v 

. 3 1 3 0 3 7 : = e * : :■ 

• 3 5 4 6 a 8 9 0 5 + C 4 

.2422478 05+ 33 

.5 05 8 P C 24E+ '4 

• 163 7 = 6535+ ’" = 

.216519685+06 




•114997615+ 71 

.=528=7845+' ' 

.25006742 *“♦•71 

# 7 

. 74 13674 E- r 4 

• 6 2 1 8 8 2 8 = E + C 4 

•1 09752785+01 


12 06 

.759037595-01 

• 28 2 078 55 5* 3° 

.30=2511=5+ : 0 

. 56 C4 : 297 5 + 74 

.26 6'’ 3 0 325+ 33 

• 5 0 5 2 9 6 3 8 E + 0 4 

.16383=415+0= 

• 21341 J2 15 + 06 

9 



• 3 1 49 9 1 65 E+ C 1 

• 95342514E + '• r 

. 271755315+91 

n 

.74 ft 20C86 = - r .4 

.620 35 8645+ 34 

. Ill 032965+C1 


12 0 7 

.772159135-01 

•29238=735+ 1C 

.3 052 043 :e + 0--. 

. 56 64 3 661 E + 34 

.2 = 2421225+ '3 

.5 04664 = ^5+ 44 

•1638865 7E+ 09 

•21010994E+C6 

= 



. 11 499580E+ 21 

•953935455+ ”5 

•29552541E+01 

n 

• 7 4=298 755 - ‘i 4 

• 61? 77C34E + C4 

.112389695+01 


12 OB 

•75536117E-01 

• 3j238139E+u-J 

.333=C684e: + :? 

.57282 '18=+74 

• 3227 0995 5+ 33 

.59401044 = + 04 

.163938675+09 

• 2 Q662485E + 06 

Q 



. 11498 .,325 + j 1 

.9544 3041E+ :: 

•3224466 35 ♦01 

. 0 

.740451 j 9 5 “ 4 

.61710694=+, '4 

.113832665+01 


12 39 

.74 3 3 935 5 E- 31 

•312 ==5245+ ~ 3 

.296330595+ 0 

. 57953 “6 7 5 +04 

.356367175+03 

.5 33319 07E+ 34 

• 163 99545E+ ' = 

.202913615+06 

3 



•114974915+ 31 

• 9 5 4= 11 5 1E+ • (] 

.351952:t5«-71 

m ~ 

.740693225- 34 

.615356535+' 4 

.115369215+01 


121 1 

• 724 5 7 34 8 E - 01 

• 32 15293 0 E+ 0 0 

.29 13925 *E+ r 7 

. 586766375 + 04 

•3=5313595+03 

.502533=35+ ’4 

. 164 058 04E + C? 

.193979 065+06 

9 



. 114=6=635+ Cl 

• 5553 8 Cl . 5 + ".0 

. T 854 281 A 5 + 01 

r 

.741 C8798E- '4 

•613E0347E+G4 

.117014585+01 


1211 

.709520025- Cl 

•33373922E+ Of 

• 286- 684 -E-*- . 0 

.594399655+74 

.439436995+03 

•5018053CE+04 

.164 12949E+C9 

.194747375*06 

9 



.114?644=E+ 01 

♦ 9558 37425+00 

. 422432465+ Cl 

.0 

.741755 15E-C4 

•61154872E+C4 

.118774925+01 


1212 

•694813535-01 

■ 33973S65E+ 3 0 

.28 3388585+ = 0 

.602515715+04 

.4= 3=91615+ 03 

• 5 0 C 9 6 4 2 5 E + " 4 

.164Z0451E+ f 9 

.190249715+06 

9 



•11495947E+ ”1 

. 9 5 6 2 3 4 c, 3 E + 6 0 

.465875375+71 

.4 

•74290C95E-04 

•6 09 446675+ j 4 

• 1 2 6669885 + 0 1 


1213 

.630415235-01 

.34654 9 46 = * v 3 

.2741 391 ? E * 0 0 

•5111438 55+ 04 

.5495 3994 E+ 03 

• 5 C 0^43 D 65+ ^ 

.164 2769 4E+ r, 9 

. 185334375 + 06 

9 



•11 495455E+ 01 

.9567 227=5+ -v 

.5137933^5+01 

. 0 

.744 = 85 52E- ^4 

.6 3715634 = + : 4 

.1227111 05+01 


1214 

• 66629 2 32 E - 01 

.35719 199E + 0 0 

•267461=85+00 

.620282915*04 

.6198=1495+03 

.499^3^ e nr+ 

.1643471 0E+C9 

.180398775+06 

9 


* 

. 11494972=+ 01 

.957152895+ 7 0 

.5707 ’1345+01 

.0 

.7 48 = 34 7 OE” ’4 

• 6 v'i 4 6 5 7 6 05 + 04 

.124915685+01 


1215 

.6524127 CE-01 

.365635255+ 0 f 

.259984755+ ~ 7 

.629858275+04 

•6=9211 00E+ 33 

.4 = 791 366E+ '4 

.164 3 = 7515+ 09 

. 17428325E+66 




.11494493E+01 

•957575875+ . 0 

.633451725+01 

.0 

.7568 = 46 2E-' , 4 

.601886365+04 

•12727733E+C1 


1216 

.63874812E-C1 

• 374^471 CE + 0 0 

.252.J04 17E + ' - 0 

.6397224=5+04 

.756068265+03 

.4 = 6655 76E+ 04 

.16442D27E+=9 

. 168 28976E+C6 

3 



. 11494 :!31E*L 1 

.9579525=5+ ' 

. 7C935774E ♦ C 1 

c. 

.774 199 16E- ' 4 

.5988071 15+'' 4 

.129799535+01 


1217 

.62527 72E-01 

.36 229 441 5 + ? 0 

.242=659? = + 0 

.64934 = 135 + r 4 

.= ■'51 = 45 35+03 

.495248225+04 

.164 373,295+ 09 

•16113C335+C6 

9 



.11 493571 5 + 71 

• 9 5 3 4 0 3 8 = 5 + J ' 

. 79 T 59 0685 + 01 


.8 153978 OE- '4 

• 59 5 36 = 6 -35 + 0 4 

.1323837 15+01 


12 18 

•611954295-ul 

. 39 -;4 4 3 21 5 + 7 ii 

.233 3 034 ’ 5 + 

. 65829762 = + 04 

.1 C378555E+ 34 

.493638615+ 04 

. 1642C588E+09 

.1537433=5+06 

9 



.11493116 E + 01 

♦ 9588 10-555+' T 

•595956385+01 

. 3 

.926650975-84 

.591458545+04 

.135003365+01 


1213 

.5987738 15-01 

•39 8 5 03 82 E+ C 0 

.22252071 E+ 0 0 

•5644=9 ORE+^A 

.1 1884279E+ 04 

.4918 = 1655+ ->4 

.163828435+09 

. 145607415+06 

9 



• 11 4926665+ 01 

• 55921333E + 00 

.1313989=5+02 

r 

•12314367=-= 3 

•597233145+04 

•13723402 5+C1 


12 2 u 

.5957 0521E-31 

• 4 .'1 6 5 0 S 9 7 E + 7 0 

.210548=75+' C 

.5670=2245+04 

.135512935+04 

.48 = 828 "4 E+ 3 4 

.163 ‘7=415+09 

.1366253=5+06 

9 



•114922195* 01 

•959612865+ " 

.114827=35+02 

. * 

.271925915- '3 

.582270555+ "4 

.138973635+01 


1221 

.572725195-01 

.414448=15+ 90 

• 1986 G2965+ 7- 

. 55 77 2 297 = ♦ 04 

• 15543964E+ 3 4 

,488 3=8 = 25 + =4 

.161758345+'= 

.127923115+06 

9 



•114917765+ J1 

. = 6 : 1 C 1 4 E '• n 

.132967425*02 

.0 

• 1 12342 94E- 2 

• 577 = 08615 + '‘4 

,138489475+01 


1222 

.5558111 OE- 01 

.42235 151 5 + 00 

.18276*06 IE* 

•63196 Q 345+04 

.1 6712759E+ 14 

.485715 01E+ -’4 

• 159 "'2 oc '45 + '9 

.116619215+06 

=? 



.114913355+ 1 

.96 04 C5575+ . 7 

.14813 08 05 + 92 

• ' 

,956875665- 32 

• 572443=65 + '-4 

.134583825+01 


1223 

.54654 >7 0E-G1 

.43 J22737E+ *• 0 

•174366325+ . 0 

• 525 = = = 7 05 * r ' 4 

.1 74721=05+ ‘’4 

.484840175+ ^4 

.152 573475+: 9 

♦11 0866635+06 

3 



•1149D695E+ •. 1 

. = 6 ' 7 9 9 = 1 5 + • 0 

. 18 37&Q°75*-‘’2 


.547456 76E-'31 

.57=33=56E+04 

.114317955+01 


1224 

.534 992375-01 

.43608991 5 + j5 

. 1 4 e j 4 : 4 5 e + : : 


r 

.4=7 224 755+ C 4 

.144244595+"= 

.975658545+05 

2 



.1149 0 455 E + '*1 

.561193355+ r 

• - 

• 

• 7 4 0 .32 fj 0 = - 4 

• 5 9 = 3 2 0 0 5 + ' 4 

• 


12 25 

. 1 5 6 2 5 3 j E + j 0 

• J 

.31 ‘724165+ 

.54758 1315+04 

. 7 

• 5 3 4317965+ ^4 

• 16346? 095+ 9 

.221309645+06 

1 



.11543376E+ .U 

•922110 6 85+ ' r 

• -* 


.7 4 0 030 455- ’ 4 

.643 .51665+' 4 

. 1 C 64674=5+01 


1237 

• 1 1 5 C 7 2 55+ CO 

• 22S?::181E* 10 

. 3 :'=9 42G6£ + 

.554=47355+04 

• 2 8^8679BC+ ;3 

.5 085 38 4 = = +- 04 

. 16363 = 7 3 =-+.= 

.216348235+06 

9 



• 115151 77 E* *1 

.941455375+ : 0 

.09511 7=35+01 

* 

.74 '033 285- '4 

.62° 759 ’85+ .4 

• 1 '9192 0 = 5 + 01 


1249 

.6453 , 646E-01 

•35238=365+ /• 0 

.255524275+ . ■ 

.535842975*04 

.7= 46 T 645 E+ 33 

.4 = «4 21 ‘4 5+ '4 

.164426235 + ’■ 9 

. 171546135+06 

9 



• 1 1 5 1 ' 6 7 E <■ v 1 

.951 7 24 ,Et ■, 

. 71235 33 = = + = 1 

r 

.7 4= 366 3 15- : 4 

.6 3 = 1= 375 + ; 4 

.125563815+01 


1261 

• 1-75 . O j G E + 0 3 

• . 

.29 :?5c:ce+ 

• 5 74 8 2 = 695 + ”4 

• 

.515210435+ "4 

• 163 451535 + ' = 

.2 7745 = 725+ 0 6 

1 



.1156245 6 E+ 01 

.? 09 4 21775+ * . 

• . 

. 0 

.74 0300 465- 4 

.64523154 = + 14 

.111571835+01 


12 73 

.1 51595255* vO 

•22989.75=+ -t 

.293515 

.57637873 5+' '•* 

.4 :314C13E+ ^ 

.5 ‘.8 8 55 3 35+ 3 4 

• 16. 36 5 4 = 25 + * 9 

. 2 ’4=^3135+06 




•11528321E+ "■ 1 

•93234848 E+ ' 

# % « 2' 1 = 5Pl r + , 'l 

e 

.74 : 034695- • 4 

• 62?5’-631E+'4 

.1 135463=5+01 




9T-V 


ITERATION NUMBER 2 0 0 0 


D T I m l 


2*76 7 2 " 21-'. S 


MODE x y 

12 ? 5 • l' l Q677245+CG . 35 7 j3 1 3 3 5 + JC 

•115094*15+6 i 
1297 .21375*005+00 .9 

.115515415+01 
1319 ,152556775+CU . 23 33 76 66 5 + G j 

.115413 >5E ♦ 3 1 
1321 .116417765+00 . 362*45 125 + 6 P 

•11519 0535+ 01 

13 33 .25 00: 0J35 + G0 . 

. 116006235+ .1 
1345 .173321595+03 . 23 733065 E + 0 ■* 

• 11554 1 iifiE+ Cl 

1357 .131691295+03 .368595455+30 

. 11 52? 3 79 5 + Cl 
1369 .232417415+00 .C 

.116204125+91 

13 31 .2 167 022 55 + 00 .25 0156135+03 

•115741845+01 
1393 .164333675+30 .388507635+03 

• 11548617E+ J1 

14 05 .3143 1451E+C0 .0 

.116402315+: 1 
1417 .239582905+00 • 26 28 8 F 71 E + C 0 

.115942 £,2 5+0 1 
1429 .197976 365+20 .4J827P73E+C0 

. 11568S55E+ 01 
1441 .34722222E+0C .0 

• 11659991 E+ 0 1 
1453 .27246256E+ OC . ?7552 G 325 + u 0 

. 116143*05+ 01 
1465 .231115445+00 . 4 2 791 J 21 E+ 0 0 

•1158 9 C 935+ ">1 
14 77 .37962 = 535+00 .0 


1489 

.3 ^5 344215+ CO 

.116797605+01 
. 28 8 08 ;95 5+ 9 0 

15 01 

.264260925+00 

.116344155+ U 
.447403495+ 0 0 

1513 

.412 C3 7 335+30 

. 116 093 315+ C 1 

• i. 

15 25 

.333224375+00 

.116986 73 E+ 31 
•3 J 0544515+ 0 J 

15 37 

.2974 332 15+ OS 

.116544975+ . 1 
.46676 •: 015+ v * 

15 49 

.444444445 + JO 

.116295695+ '1 

« iJ 

1561 

.3711 C 5525 + OC 

•117157925+ 61 
.312 = 2, 92 5+ 00 

15 73 

• 33 J545595+C3 

. 1167457=5+ J 1 

1535 

.4766 5 1 35 5+ '0 

. 11649? .,75 + : 1 



.117323875+ =1 


PHD 

cv 

•235381615+ t r 
. 94 5 4 79 8 -''5 + - ' 
.27 ■ U 9 4 C 3 5 + •• r 
.59 7 ; 4 3 4 7 £ + : 0 

• 27569 C. u,2E+T 0 
.92 35 068 6E+ . 

• 2 1 3 2 1 7 6 4 E + ■ ' 
•93875322E+:: 
.24939397E+ jC 

• 86496465E + ',0 
•258107115+ CC 

• 914°3£,?3e* "5 
.2330122^+0 9 
.9323 j 8 2 1 5 ♦ 0 0 
.23 01 1954E+ :.o 

• 87274346E* 9 0 

• 22° 7 9366E + 9 9 

• 9 317 79 585+ r ’0 

• 1586 215 4E + -JC 
.91359373E* 00 
•2119£41 C E+ r ? 
•860821675+ G 0 
•203235455+ C C 
•888?5S4°E+( C 
•123011635+ 10 
•9652331 'E+ JC 
•194737165+ iC 
.3*9188575 + •? 3 
•1793923 9E+ 00 

• 876461 4 E ♦ OC 
.1 '6013 715+ 1C 
.39222758E+0C 

• 17842549E+ r 0 

• 337 Q 3392E + C f i 

• 1584 1891E+09 
. 3 6 4 2 7 5 3 2 E + i- 0 
.933497165-1-1 
. 3 7954914 E+.-f 
.163 39 "0 05 + 00 
.3275 0759E+ DC 
. 14 02 054 9E+ 00 
. 8 5 2 3 8 9 9 5 E + 0 C 
. 7 9 1 S 2 5 3 9 E - 01 

• 8 6 7 1 9 0 3 6 E + ■; * 

• 14 Q 936^4E + ' C 

• 91391393E + -IC 
. 1244 573 = E+ •. <1 
. 8 4 J 7 9 4 1 ; E * J C 
. 7 1 2 1 6 2 4 l E - 1 
. 8 5 5 1 3 9 4 3 E + ' f 
. 1 3 6 2 4 1 5 ~ ♦ r - 
. R 1 '48 qqc E + ; 


'J 

T4 XY 

.6596-5. ' 1E + C^ 
. 95962237E +C1 
. 6344 1 34 3E + C4 

• ^ 

• 59S76 764E + 04 
. 52352*385+01 
. 67956831E+C4 
. 116234695+02 
. 6347? 644E + r 4 
. 0 

• $2*43 2995 +84 
. 626 07 436E + 91 
.7.3229312^ + 04 

• 14 L'l 9 9 06E+ r 2 
. 664 7641fcE + 04 
. C 

. 6627677=15 + 04 
.845219145+01 
.74378601E+04 
.1961 C 5125+02 

• 69297217E + C4 

. 0 

•699079 03E+04 

. 1 C667212E + 02 
. 76876 427E + 04 
•245765315+02 
•721740365+04 

• M 

.73221 165E+04 
.127693855+92 
. 784 08 829 E+ 04 
.2824 1 111 E + 02 
.750468975+04 

• ti 

. 762065675+04 

• 146925865 ♦ 02 
.79471 0155 + 94 
. 3 0 89 0 3775 ♦ 02 
•779 0 39 8 PE+ 04 

. 78 88 2 1 .’75 + 04 
.264177705+32 
. SC? 7 ?9755 + 04 
•326591965+02 

. 9 C6°6 66 1 5 + C* 

.8 12946195+- 0 4 
.1794 O 1735+32 
.512671315+3^ 

. 3 368 76 = 5 = *92 
.6337156''=' + '"' * 


T«r = 


346735° =5-02 


V 

T H X 7 

.11153433 5+0 4 

• • 

• - 
* 

•548631425+03 

m 

•13973545 E+ 04 

9 

• -■ 

.685050295+03 

• 1753 6 052 E+ 04 

. r 

n 

• ■ 

.9® 486 003 E+ 33 
• <-. 

.267933595+ 04 

• 'J 
* 

9 V 

• z 

.1 3 1678065+ 3 4 
.35 158 69 4 5+ 04 

r, 

• - 

r\ 

.0 

.16594302 5+ 04 

.42 1149095+ 04 

• - 
• ■ 

*1 998187SE* 14 

r 

.47544294E> 14 

• - 

• 23242853E:* 14 

• * 

• 5151833 IE* ?4 

• - 

r 

• . 

.26 32 -4 37 5 + ^4 

O Z 

.541865595+04 


SOS 

“U 

• 497" ; 37545 + ‘ 4 
.74 5 133 o.fe e- 04 

• 5 I 57D.069E+ °4 
.74 f. DOC 485-94 
•508 8 43385+ 04 
.74 ''036685-34 
.496 024 515+ 04 
.745363 3 4 E - 3 4 
•515944195+ 04 
.743000 01E- fl 4 

• 5C861051E+ 04 
.740000875- 04 
.494220475+ 04 
•740124 C 1E- '4 
•5161 73 175+ °4 
.74 0004 015- 04 
•507665265+ 04 
.740000 Q 9E- 0 4 

• 4 89 4 04 8 25+ 04 
.74 0126 085-04 
.516229625+ 04 
•74C3C0015- 04 
.506390635+04 

• 74 n 0 01 035- C 4 

• 48 52306 0E+ 04 
.740119° 45“ "’4 
.516097235+ -4 
.74 0000 025-34 
.5 04947 ‘■’55+ 04 
.74 H 0 0 1 04 5- n * 
.481868355+04 
.74nil294E-' , 4 
.515745295+04 
.740000 025- •'4 
.533442985+04 
.740001 " 1 E - C 4 
.479382655+ n 4 

• 7 4 .0-1 06 5 1 E- 0 4 
.514978955+ A 4 
.74 OiJOG 02 5 - ij4 
.5019 09 275+ 04 
.740003985- 04 
.4 778 05885+ 04 
*74 0101 65E- '*4 
•513883835+ 04 
.74 0003 015- 04 

• 53 0*41 '45+ 1 4 
.74 09 909 3 E- '. 4 
.477361335+04 
. 74 ;! '"SB 7 8G- 3 4 
.5 12725 155+ ’4 
.74 ’’330 :n r -0‘4 


r 

T 

.16456463E+- Q 

• 6 i '06562*5+ n 4 

• 16346 0 2 9 E ♦ n 9 
.64 6 ° 06 11 5+ ' 4 

• 1 63 63 '6 ;E + «. 9 
.62°44;77E+ 04 
.164 67 78 SE+r-9 

• 59° 1 8 o 9 E + 14 
.16351033E+C9 

• 646 °629 0 E + 0 4 
*163793925+39 
.628829435+04 
.1648235 ‘ E + J 9 
•5 q 3?17215*04 
•16354494E+C9 
.6474 81 5 C3 5 + n 4 
•163 a 0022E+89 
.62683915E+.'4 

• 165 17553E + 0 9 
.5822*9525+04 
•1635921 15+09 
.6475677 lE+u4 

• 164 ';31 = 35+c9 
.6232423 3E+C4 
.1654554 i E+ ;.9 
.572 3 0 98 95+ 04 
.163653365+:° 

• 647188335 + '. 4 
.164 172 3 2 5 + 09 
•61964C36E+04 
•16569146E+09 

• 5643552 OE + 94 

• 16373641E+'*9 
.645242955+ 64 
.164322085+09 
.615930335+04 

• 1 65 8824 5 E+ 09 
.55° 4 Q 8 6 j E ♦ ' 4 
.163832925+99 
.644 .'6*5 7E+ 14 
.164457565 + <■■ 9 
•6121 03 335+ 0 4 
.166 0 3 7 1 6 E + v 9 
.5547819 3E+ 04 
*163929335 + -j9 

• 64 '661275 + ‘4 
.164611175+ =9 
.6 3472315+ .4 
.166131115+ “9 
.553' *i 5 9 . 5 + *" 4 
•164 ' 528 .5+=° 

• 63711345E + .'4 


P 

M 

•1573332 5 E +06 
.13*606975+01 
. 1927698 :.E + 06 
•117203345+^1 
. 192233625+06 
.118164*35+81 
. 144668615+06 
. 139S713 4E + C 1 
•1778 7 J° 55+06 
.122033935+01 
. 179609475+06 
. 123508955+C1 
.131711735+06 
. 146464115+01 
.163990155+06 
.128651435+01 
. 159 03675E + 06 
.131985645+01 

• 1 0225007E + 06 
.161537695+01 
. 15;)82828E + C6 
.1342372 :E + C1 
. 1397 04995+06 
.1*0478975+01 
.309746865+05 
.174215565+01 
. 1282638 5E + C6 
.139845755+01 
. 1224 T41 3E + 06 

• 148684955 + 01 
.66 C17049E + 05 
*1347 04865 + 01 
•12629546E+C6 
.145511555+01 
.107264345+06 
.156487865+01 
. 555877575+05 
.193180225+01 
.114911765+06 
•151265775+01 
• Q 4192622E+05 
. 163844625+ul 
.483129965+05 
. 199804785+01 
. 104341245+06 
.157332895+01 
. 5 298 133 • 5+05 

• 17074 7° 45 + 01 
.432417925+05 
.2-'477 a 236+Jl 
„°472733 35 + 35 
.1626 0 47 8 5+ 01 


IB 

9 

9 

9 

1 

9 

9 

1 

9 

9 

1 

9 

9 

1 

9 

1 

9 

a 

1 
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U 
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.4 3986195+07 

. 325211 : 35 + 1 : 

•11 69 442? E+ =1 

. 11 :?i®2?5+ 
.429669465+' 1 

.335158 5 4 5 ♦ Y 4 

. 1? ?6 4 = 72^ ♦ 02 

16 09 

•363657975+00 

.515 9682=5+ :■ 0 
. 11670 „45E+ -1 

.65453 c 655 - 1 
.543365265+0 

•®2268537 r ♦04 
. 34 1 43 4 "55 + 02 

1621 

• 5 19259265 + 00 

• 

. 11749991 E+ Jl 

• 124293 2 -"5+00 
.8022311 15+ = 

• 35930722= + 04 

r- 

• it 

16 33 

.436866335+10 

.33741574E+ : 7 
.117117635+ vl 

.99 " 496955- 1 
. a 2 ' 9 0 7 ° 4 E + - r •" 

. 5556123°5+C4 
,2" 4 -'91945 + 02 

1645 

•336339365+00 

.524 12 2 78 5+ ' 1 
.116972935+11 

.61 12 72935- " 1 
• S3 1? 171 3 5 + 1 ’ 

. 53372562=+ J4 

.34172°635 + J>2 

1657 

.541666665+ Du 

• 

.11767 '965+ <U 

.113651445+-:* 
• 79 4 1 324 J E ♦ ' 0 

• 38377 851= + 04 

• ■ 

1669 

.4&974748E+ u0 

. 34 95 3 569 5+ 1 . 
• 11 7291375 + >Jl 

.?83972365-=l 

.612322385+00 

.3744 J8425 +Q4 
. 213=38235+0? 

16 8 1 

.429972745+00 

• 54 2945°5 5+ O 1 
.1177? 144E+11 

.5766121 ''5-71 
.822731125+ J 1 

.84541 7375+04 
. 33=577045+ 02 

1693 

.574 07407E+ CO 

. 0 

.117842 .05+ =1 

. 104 C 0311 E + :i 
.786189245 + 0 0 

.9i0715 0285 + 54 

. P 

17 35 

.5 '2629145+C3 

.361572315+00 

.117464915+01 

.SCI J8285E- 01 
.6 039 07° 4=+ •*:• 

.391841=55+14 
• 2223664 3* 5 + 02 

1717 

.463115125+ JC 

.561539095+ :P 
.117256375+0 1 

• 54782653E- '1 

. 3 1 4 : 4 r 2 3 5 ♦ 1 1 

.357447795+14 
i jjcc'-ooar + T2 

172? 

.6 >6461485+00 

. 0 

.113013045+01 

.952716145-01 
.77°3?72 15+ cr 

.52=488315+04 

. 0 

1741 

• 5355 u87 95+ 03 

.373525855+0= 

.117635455+01 

.724737955-11 

•795659165+01 

. 908141765+14 
•22=5467 15+ 02 

1753 

.496257515+00 

• 5311035 0 5+ ■'’7 
. 11743 1295+ 1 1 

.521717115-01 
• S 1 5 5 2 5 7 5 5 + 10 

.36=154065+04 
. 33 2 09:144 = + 0? 

1765 

.6338838 9E+ 00 

• ^ 

.11? 13 4 09 5+ U 

.87356746 F- Jl 
.77*752125+17 

. =5 18 1 0865 + 04 

.0 

1766 

.63033 :935+ 00 

.467835665-01 

.11813®925+01 

.889179155-11 

.772756665+07 

.940953585 + 04 
.247242975+01 

1767 

• 622498095 + 00 

.896=31 71E-01 
.113097535+ Jl 

.5911 ?44 = 5- " 1 
.774600655+ . 0 

• 93 741 3? 95 + 14 
. 4 Q 393?97=+01 

1768 

.615290275+00 

.129105= jE+°G 
• 11805?54E+’Jl 

• c £4? 36£ =5- Jl 

.776 3 0477F. -1 9 

.9 33 833645+04 
• 7283 1l4 7 5+''l 

1769 

.6.3624585+10 

.1654446 :•£+>??, 
.118024565+ 01 

.3 74. .-319 15-01 
.777887115+00 

•932777335+04 

. 936=7 676= + 01 

1770 

.6:2432255+00 

.19929647E+ J 7 
.117991675+ 01 

.356 1433 15- Jl 

.779362665+ 10 

.931216045+04 
.11354= 07=+02 

1771 

.596654895+00 

• 23. <679395 + .« 9 
.117961 135+ 61 

.6315662 = 5-91 
.78 17442 "5+ 0 *; 

.93 !24r>445 + 04 
• 1333 "8835+02 

1772 

.591243375+10 

.26 462425+ ID 
, 11793262E+ u 1 

• 0 16 362195- Jl 
. 76 2 ‘ 425 - E + ) p 

. 929=486 8 5+04 
.151229385 + -02 

1773 

.536155655+ uJ 

. 26 827531 5+20 
.1179 05775+ 11 

. 7 7 ° 4 ;)7 2 6 5 - 1 1 
. 7632669 ‘ 5 ♦ 0 

.929=41155+ 04 
♦ 1 6? 2 3 065 5 + 02 

17 74 

.53 1 3556 45+ DO 

.3145154! 5+ 1 ;■ 
. 11768 *435+ '1 

.747797695- >1 
.7544 25 455+1" 

. 9282347 "5 +14 
• 1 8 6 7 7 2 4 8 5 ♦' 2 

17 75 

.576312125 + C J 

.339353335+ J J 
. 117656455+ 51 

.719 171145- n 

. 76 =5 2 K 1 25 ♦ ■ J 

.9 27 05 84 3= ♦ J4 
• 20 32 =, '24 = +"2 

17 76 

.572497965+ Du 

.362937*0 = + 01 
.117833635+ 01 

.£37369695- -■ 1 
.736572015 + 1" 

.925411325+04 
.22 0 023525* 02 


345725= ’ 5-12 


y 

tpvz 

2° 1QOP77 = + -»4 

5 5 75 2925 E + *4 

r 

2 1 “35551E+ 74 
55 5 1 2 45 1 E * 74 

34256825E+ J4 

i 

56=33344= * 24 

>n 

36461912E+ 04 
n 

5? 949281 5 + 14 

7 
i *> 

3?^5 3565 E + 14 
*» 

54 97 5 944 E + 74 
4 ''629 27 IE* 03 

r 

61 84 6 575 E+ ' 3 

i 

11942894=+ 04 
1=3914655+ 04 

T 

1 8 7 0 0353 5 + 94 

2 2 ^4 29 <02 = + 0 4 
25131972E + 74 

7 

282642435+ 14 

r 

31377 93 4E+ 74 

3 ° 3*5 33 1F+ .14 
373934645+ 74 


SOS 

MIJ 

,49?fil?33L + 74 
.74^010395- '4 
.4 7 71 05 05 = + 74 
.74 0 796 5 IE- 0 4 

• 5 1 1 51 7-665+ 04 
.7 4 7 1 07 - iE- 74 
.4=6969°? 5+ 74 
,74??'1P845- 74 
.4775614 5l+ 04 
.74 7 795 62E- 74 
.51 =269 71 E* "4 
.74 ^ 900 7 1 E - v 4 
.4=>522345E+ 14 
.74 0 310 8 15- 74 

.470 71373E+ -"4 

• 7 4 0 795 26E- 14 
.5 08=85315+ C4 

• 7 4 ° 0 13 1 1 E - ** 4 
.403575 135 + 14 
.74-901=765- 74 
.479638 19E + 74 
.74 CQ9554E-'4 
.5 07676 195 + -04 
.74 rone hie - 74 
.492022° 0E + >7 4 
.74 : 2:3 72E- =4 

• 4793 96 285 + ""'4 
.74 ??95’??5- 04 

• 5C634=29E+ f '4 
.74 orr<0 'lE-*4 

.507472115+74 

.740103 92E-04 
.5 16525 345+ 0 4 
.74 1013 03 E- ?4 

• 5 7 56 76 7 7E+ n 4 
.742170 14 = - '4 
.5 14596 78E+ ">4 
.74 0 100 O5E-«0 4 
•50332622E+14 
.74 7 010 1£E- '4 
.501 7449 3E. 14 
•74 9 33311 E- *4 
.50 717963E+' 4 
.74 " 919155- i4 
.4QPJCI! "6E+ '4 
.74 -2-’227E- '*4 
»45t*oa 1 0E+ ; 4 
,74 *9 00 275- -74 
,4=46:1 a SE+ 14 
.74 ’973775- 4 
. 442560435+ 0 4 
.74 7 7 00575- 14 


E 

T 

. 164 7 D357E+ ' = 
.6 ''443276- + 1 - 4 
.166 1 Q 95 'E + o? 
.553 '59 0 1E+ ^4 
.16416754E+19 
.6334591 5 E + 14 
•1648647=5+ ’9 
.599 32 74 3 E+ >'4 
.166 l9223E + '-9 
.554 169485+14 

• 1642714 r E+'‘9 
.629 71595E+':4 
.165 019 6 jE + 79 
.594 495 7 5 E + .'4 

• 16619129E+.9 
•554573245+14 
.164 362=25+ 09 
. 6259073 CE+04 
.165 !65?35 + -"9 
.589923® 35+ 14 
.16619 84 r E+ 05 
•55546744E+14 
.16444 03 75+19 
.622 15 C 22 E+ 94 
•16539373E+1? 

• 58567391E + C 4 
.166 17683E+-‘9 
.55642 1765+ 04 
.16450465E+ 09 
.618 16645E+ 4 
.164528935+09 
.621 179665+04 
.164316305+19 
.615=1 322E+04 

• 164 2 1556E+ ? = 

• 6 1 6 ° 1 £ 4 6 5 + . 4 
.164377845+"° 
.6144837 1E+ n 4 
.164443995+^9 
.611514645+04 
.16451744E+:9 

• 6 177 9 9665+04 
.164695 715+ 09 
.6=4 1051 1E + ‘0 4 
.164 795 79 5 + ' ! ° 

• 5 9 9 ?, 8 4 7 : 5 + ; 4 
.164=53435+ r ? 
.595436955+74 
.165112175+09 
.59098 C375+ . 4 
.165264565+79 
.586192625+04 


P 

W 

.73284*625+05 
.1773584=5+01 
.3968 128 4E+C5 
.208343865+01 
.°6 u25132 = + J5 
. 1679917 1E+01 
*6492798 1 E + 0 5 
. 1836°7£6E+: 1 
.370649665+05 
.211 G 11875 + C'l 
. 781622945+05 
.173198895+01 
.577722865+05 
. 189635255+ Cl 

• 3 4 9 7463 6 E + Q 5 
.213233565+01 
. 71 ,-1643295 + 05 
.178229175+01 

• 5163816 jE + 05 
.195217885+01 

• 3 32 6731 2 5+ 15 
.2150558 15+31 
.64669S7°5+G5 
. 193086345+71 
.46355:745+05 
.203445615+31 
.317229675+05 
.21671 C8€E* Cl 
.58909115-5 + 05 
. 197777665+01 
.672442475+05 
.185592535+01 
.601706795+05 
.185771595+01 

• 5°55699 75 + 05 
.196207975+01 
.586 0554 oE + 05 
.187355835+31 
.571266095+05 
. 1 88 7C.6 375 + 01 
.551592715+05 
.190535125+01 
.531669115+05 
. l a 25931°E+ul 

• 5 = q 6®52 2E+ 05 
. 1 9'4 8 2 5 7 7 5 + 0 1 
.486042535+05 
.197346955+01 
.463? 175 1E + 05 
. 199685575+01 
.439861 715+05 
.202636015+01 


I? 

a 

1 

9 

5 

1 

Q 

9 

1 

9 

9 

1 

a 

1 

a 

9 

9 

9 
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5 
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ITERATION M U M 8 3 3 2 CG ' 
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5 - .2476 72 A 2 5 
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y 

Y 

RHO 

U 
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CV 

THXY 

1777 

.563389455+02 

. 3353 973 C 5+ 0 j 

.65514 = O'. 5-, 1 

. 3 234 8 0425 + -'4 



• 117S12cOE+-U 

.78 75715 ’5 + 

. 235629955 ♦ 0? 

17 7S 

•564465735+ OJ 

. 4, )t 846935 + o ) 

.6261 06° = 5-;.l 

. 92 08= 8965 +u4 



•117791 295+01 

.783526315+ *, '• 

•251116355+02 

1779 

• 5 6 * 7 0 8 4 u 5 + 00 

• 42733 7 : j 75+ ' C 

• 6 '.*' * C 7 -* 2 7 5 - . 1 

•91786643=+ "4 



.117771465+. ’1 

.769446655+00 

.265724°* 5 + ^2 

17 3 0 

.5571 'll 25 + jO 

. 4471 16915+ CO 

.573272535- 01 

• 91429 0345 + 04 



.117752425+ ">1 

.79 03 3C245+ O' 

.279648265+0 2 

1781 

.5536293 15 + 00 

. 4 6 6 8 6 1 5 5 + 9 O’ 

• 5 4 9 0 2 2 4 0 5 - i 

.31 0274 01 6 +'■'4 



.117734 095+ 01 

.791182445+ -r; 

• 29251 °5 15 + 02 

1732 

•550279925+00 

. 484 3° 6 15 5+ w 0 

.527649125-01 

. 5 65° ° 7=9 5 ♦ C 4 



.117716425+01 

.792 0 06275+ u 0 

• !04 ; '125°5+ r: 2 

1783 

.547941211+00 

.5021 51135+1' 

.51 T3 9:70 5- '1 

.3C1665485+C4 



.117699325+91 

.79 28 C44 4E+ 1 0 

.313667665+02 

178a 

.543932555+ 90 

•518259145+00 

.4978 75585-01 

• 9 9 75 7 9® 75 + ft 4 



.117652755+01 

•79357=435+ K 

.321 ' , 14°75 + 02 

1755 

.5408543 GE + 00 

• 535922915+0 0 

.49 C-9 32 8 0 5-01 

.393932935+04 



.117666675+ 01 

•7943334 9E+ O' 

.325692525+32 

1786 

.537367675+90 

.55214 05*5+ 'JO 

.469 J16 09 5-01 

.99 '65 2 78 5 + r 4 



.117651 015+ 01 

• 73 5-0 686 = 5 + 'G 

.727354495+02 

1787 

•534994535+00 

.567956125+ 0 0 

.4 91 1409 4 5- "'l 

• 99 756 4705 + 04 



.117635745+ 01 

• 79 5 786 9 1 E + C ' 

.329265275+02 

1738 

.532167615+00 

.58341 •; 235+0 *• 

.4°44696°5-'l 

.884432875+04 



.11762 0 8 2 5+ 0 1 

.79 6489925 + O' 

• 3 27° 5 7925 + 02 

1739 

.529399395+ 00 

.59 3 54 0 44 5+ 0 

.498 j 03 0 05 - 0 1 

. 33 0 9 5 2 J 2 5 + 04 



.117606215+01 

.797179365+ ?r 

.327777845+02 

1790 

•526655C35+DG 

.61333 1685+ 0 *i 

•50015646E- U 

.3767° 0215 + 04 



. 117551655+ ul 

.7973 56755+ •' 

. 3280' 3045 +02 

1791 

•524G17065+C0 

.627966535+ <: j 

.5C17C5B5E- 01 

.9715*9135+04 



.117577805+01 

.79=523515 + *■* 

. 328774725 + 02 

1792 

.521 3 9 <-3 5 E + 9 0 

•642325925+ 9 0 

•5C1451Cl5- n l 

• 8 6 43226 = 5 + 04 



. 11756 394 5+ 01 

• 7 9 9 1 3 1 0 ' 5 ♦ 9 

. 329964675+02 

1793 

.513799675+00 

.b56483375+ 0 0 

.5C1485795- 01 

•3537*2505+04 



« 11755 0265 + 91 

.799836515+ : 0 

.331331215+02 

1 7 9 A 

.516239915+00 

• 67043 175E+ 2 9 

.499015355-01 

.937128415+04 



.117536755+01 

.6004 73255+ 90 

. 33277 l?O r +02 

1795 

.5137^6235+00 

.664332.285* 0 J 

.497626 .46- ;l 

.81 021 133E+J4 



.117523395+01 

.60111*415+00 

. 534825145 + 02 

1786 

.511194155+00 

.69° 06 5 25 5+ 0 i 

.49 0695755-01 

. 766°?5675 + ''4 



. 11751 '12 5* Cl 

.801743135 + -JC 

. 336722405+02 

1797 

•5-9699 055+ 90 

• 7117 0 5125 + JO 

.4891 02 8° 5- '1 

. 6974 1 8155 +0-4 



.117496955+ Cl 

• 30237251E + jf 

•34=116365+32 

1793 
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. 725275735+ JO 

.482680045- -1 
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.117483855+ Jl 
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• 34 136°645 + 02 

1799 

.5j374 26 05 + 0 9 
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. 11747 1795+ 01 
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ie c: 
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.117457765+ ul 
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• 
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. 971166795 + C4 



.118355135+ Jl 
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.6 .12701 05+ .0 

.397187435+- r 

.595387275- '1 
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*■ 
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9 
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» A 5 O J A 2 1 8 C ♦ - 4 
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• 165752475 + 1.9 

.374840985+05 

9 

• *■ 
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• 2 1 0944255*01 
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•163909125+0° 
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9 

* •' 
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.567621095+04 

.213621745+01 
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.2161 J7815+C1 
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c 
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.559812925+04 
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.549661335+94 
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•166307485+ . 9 
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.7409 05 ’15- 04 
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.56 3 '7788 5+ 14 
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.3 n l736C 65+05 

Q 

- 

.74 .1597 9 4 5- "4 
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• 22 1 1°12 1 5*5 1 


.57 1 018245+ '4 
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.16646795E+ . 9 

• 297G9C6 -5+05 

9 

: 

.7 4 0C12 225“ *'4 
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.221614375+01 


.573666 45 5+ 04 

.479756 665+ = 4 

.16646=135+09 

.256048695+05 

9 

0 
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.220472865+01 


.5 69 9 = 6165+04 

.4 7°559 755+ ‘14 

.166311355+09 

.300492795+05 

3 

- 

,74 ''95441 5-04 

.556369455 + '-4 

.215396345+01 


,567251515+ 04 

.480 r -25 4 05+ 04 

.1661764 lE + -*9 

. 3 032669 45 + 05 

9 
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.557502 305+04 

.218276735+01 
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.48 0430245+ 04 

.165979715+09 

.3051 29465+05 
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. 
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.165685 955+09 

.3 J6559695+05 
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,C 

.74 '346825-C4 

.559356785+04 

.215839135+01 


,56 0994535+04 

.481163925+ 04 

.1652085 15+09 

•3C6927255+C5 
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r 

.740723=05- 04 

.560292975+04 

.2141483=5+01 


55722514 5+ 04 

.4 8171 3 685+ '4 

.164444195+29 

.307686135+05 

9 

- 

.741655375-04 
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.523194195+04 

.111581365+05 
• 263563595 + 01 

9 

24 91 

• 11262656E+ Ol 

• 91 1515495 + 2 0 
. 12 C 946725+ -j 1 

. 235539835- . 1 
.6971 73 •*‘65 + 9 

. 1 94966495 + 05 
•2737 s 9 645+ 02 

• 54 3^ 057? z+ 

.40 1551 435+ r 4 
• 1 4 T 6 77 575- " 3 

.166224555+^9 
.5519774 :e + -24 

. 141414635 + 05 
.245472445+01 

9 

24 13 

.124349525+11 

• ■* 

•12C44445E+. 1 

. 235 : 358 15- :l 
.603341165+7" 

. 1226 c 595~ + °5 

a 

-1 

•472672 375+14 
.74 1900 *15-14 

.165=91725+09 
.53 2 262? 85+ M 

.135509585+05 

.259550115+01 

1 

24 25 

• 11951699E+J1 

.59=423235+ = 9 
• 121276425 + -01 

. 18663" Q 25-01 
• 6891 15255+9;' 

.113539055+05 
• 24 "‘9 50425 + 02 

.5077° 5355 + 04 

0 

.460467 "'65+ ’4 
.74 1 993355- ’’4 

.167336285+09 
.521 5126-5+14 

. 1 058 69 9 E + S 5 
.265497375+01 

9 

24 37 

• 1 165673 5E+ 01 

.93 j93 2945+ C P 
. 12 ;17? 11E+ 91 

.226273535- 71 
.692541575+ JO 

.1 6579 264 5 + «75 
• 2 7 9° 9 £94 r + ^2 

.541365355+14 

.40 3424615 + v* 4 
• 74': 978 755- '4 

.16619- 2 725 + '9 
.551156845+ '4 

.1356295 C 5 ♦ G 5 
. 2468498 95 + 0 1 

9 



A-21 


ITERATION NUM95R 2SG0.. 5 T I M 5 = .247672025-06 TJ 

MODE X Y RH3 'J 

GAM CV THXY 

244=) .123676955+01 .9 . 2 2 " 7 06 55 - ' 1 .12 386* Q 25+ I '5 

• 12 25 2 ' 72 5+ J 1 .68 • 78 * 415+ 0 -1 .2 

2461 .123476585+01 .612114815+1 1 . 178 994 0 2E- C 1 . 1 1 45 4 5895 + r-5 

. 1 2 .*4 0 6 * 1 5 + J 1 .684566=25+ 0 r ‘ .23715144^+02 

24 73 .1206 Jv.53E+l-1 .95 -64 358 5 + C 7 .217196765-"! . 106634928+15 

. 12 >31254E+ 01 .68786473E+ 1C- . 2683 1 876E + 02 

2485 .132601415+01 .0 . 2 05= 98 4 75 - 0 1 .125913045+05 

. 12 :59£92E+ 01 .67°, 17444E + 7 0 .0 

2497 .127536555+^1 . 625 2 3 j7E+ 0 0 .17 0245535-81 . 115569925+05 

• 12 049 76 9 E+ ) 1 .6815 414 1 E+ 0 0 . 23 2° 0687*: + 0 2 

25 09 .12473543E + 01 . 97 '69 072 E+ 0 G .2 7-83818 Q E-'l . 1 07 4 9 8565+ 05 

.12 i44218E+01 .683*17785+ 7.0 . 26 55 4 48 OE ♦ 02 

2521 .136634:85+01 .0 . 19 27 4996E- 71 .126132275 + 05 

• 12 067*2.75+11 .6755 17155+00 . U 

2533 .1317J346E + 01 .638177645+00 • 1 6 3 J 64 9 1 E - 0 1 . 11659 244E +05 

.12:581125+01 • 6 7 £ 7 6 5 6 5 5 ♦ ; 1 . 2283 1 567E+ 72 

25 45 .123984355+01 . 99 112 ' - 39E+ 0 .200171305-^1 . 1 08 3 7 1 17E ♦ 05 

•12 '52 6 8. 4 5 + 01 .687575575+70 .262512195+02 

2557 .14 o‘7£719E+ul .0 . 1 9 02731 25- 31 .12722 489E + 05 

• 12 J76203E+91 .672801 q 95+ IT .3 

2569 .1360 J496E + V.-1 .65161 :235+0'. .156478725- 01 . 117582935 + 05 

• 12 66 6 74 E ♦ 9 1 .675= 3 0? 45+ >" . 2 2 34 2 878 5 ♦ =2 

25 8 1 .13336:165+01 . 1 J 1 1 5 5 IS E+ 0 1 . 19 19 05 14 5- 01 . 1 C9 258 1 1 E + C5 

. 12 061 4 G4E+ 01 .6776728 85+0.0 . 2598 4 3595+02 

2593 .145074155+01 .0 . 1685 137PE- >1 . 12 8 2 9 4 065 + 05 

.12 0847455+01 . 67 1 0 21 =55 + " 0 .C 

26 05 .140439945+01 .665354795+ 30 .15037661 5- 31 . 11857 156E +C5 

• 12 375511E+ 31 . 673 9 28 1 9 5+ . 0 . 2183 04975 + 02 

25 17 .137377005+ 21 . 1 3333273E+ 3 1 . 1 £ 3747 1 75- J 1 .11012 3535 + 05 

• 12 07 24 04 E+ C 1 .6747 3229 E+ * 0 .257432765 + 02 

2629 .1495:9765+01 .9, . 1 5 7* 5 1 5 99 - 31 . 1 2 Q 34 435E +05 

.12 0935845+01 .667169385+0 0 .C 

2541 .14532870E+G1 . 67 94 4 8 21 5+ : " .144668635 - "1 .11 95 49 38 5*05 

• 12 98 4 6 55 E+ 9 1 .67 0* 351 2= E+ 3 8 .213012605 + 02 

2653 .14255 0475+01 . 1 05 52 1 56 5+ U 1 . 1 761 12 0 OE-0 1 . 1 1 09 5 7 1 = 5 ♦ 05 

•120797165+01 . 6 7 16 55 76 E + 3 0 . 25477 14 "E + Cl 

2565 *15411 0432+21 .9 . 14 7 0 4966E - 3 1 .13 0381255 + 05 

.12102751E+ Cl .664236 1.35*00 . ■.■ 

2677 .149788235 + 01 . 69 393 " 4 .3 E+ : f . 1392 45 9 7E- 3 1 .12 3517145*05 

. 12094 1395+ l 1 .666991095+ . 2 f:76 3 46 D5 + ^2 

2569 .1473 = 7355 + -01 . 1 0 777 . 71 E+ . 1 . 1685579 35- . 1 . 1 1161 ^495 + 35 

. 12 389 3755+ ‘J 1 . 668 5 2*5 5 5 + : 2 . 25 228 3245 ♦ 02 

27 j 1 .15389 44 55+01 .0 .137289575-31 .13141 1425+^5 

• 121122645+ fJl .661212725+0 0 . 0 

2713 .154737435+01 . 7 3 384 4 97 5+ 0 ' .134:23225-01 .121477035+05 

.121 04 ."035+ 31 • 6 6 3 8 3 3 2 ' 5 + J .202262275 + 92 
2725 .1524334 -3 E + G1 .11 JjS7 015+3 1 . 16 .06 58555- 31 . 1 1 2 6 7 Q 55 5 + 35 

,12099419E+il • 6 6 5 2 ° 9 1 2 5 + j 0 . 25 0 1 3 663E +C2 

2737 .1633322 85+31 . ;i . 1 2 8 1 25 7 9 E - ' 1 . 1 324 4 1 1 35 + ' 5 

.121222235+01 .653:89575+39 .0 

27*? . 15939 748E+ 01 . 72 42 3 9 75 5+ 0- . 1 2 3 9 9 r 3 6 5- : 1 . 12 2 4 3 0835+35 

• 121142635+ 31 .66 05323 15+ ' C . 1 9 6 ° Q 13 Q E + 02 


3467359 '5-92 


V 

T w X Z 


P 


15 


5 "31 Q 245 E + "* 4 

5393 c 6935+ 0 4 
-*1 

4 9 7 4 96 14 E ♦ 54 
537244905+ 0 4 

49 0823435+04 

n 

534456115+04 

■T 

*8 3271235+ 04 
532511175+04 

474984655+ 34 
53 09 4 25 = 5 + '4 

466132665+ 04 
528694855+ 04 

*56922385+ 04 
526815485 + "4 

**7579365+ 0* 
.5 2 5 7652=5+04 

n 

9 

.* 3 332 093 E+ 14 


SOS 

MU 

•470548=2 E+ C.4 
.74 n r?5 3 -'If- '4 
.467917115+04 
.74318? 6 E - 0 * 
.4813125 05+ 04 
•14574237E-23 
.467399 8 35+ 04 

• 740000 'HE- 1 '* 
.466955525-+ 04 
.740000345-04 
.4 8 t994 BQ E+ ^4 
•740C77375-04 
.464713625+ q 4 
.74 ''030 01 5- 04 
.466021 095+C4 
.740000325-94 
.48 0179685+ 04 
•74C072 865-44 
.461984425+ 04 
.7 4 0 045 9 05- "4 
.*65193 ’ 3 E ♦ 04 
.742963865- 04 
.47929393E+C4 
.141^07755-03 
.4 592 08 3*5 + *'4 
.74 0 000 005- 0 4 
.464*54 485+ '>* 
.74 4030 325-04 
.4 7 0 4 1 5 115+ "4 
.7 4 *073 145- " 4 
.456383845+54 
.742200 105- n 4 
.4637P4°35+04 
.7 4 r ’O0O3 05- 54 
.477596*35+ C4 

• 7* 0 068 29 5- n * 
.4535 09 245+ n 4 
.7 4 0 900 9 05- 04 

• 463162 55E+ '4 
.740039315-04 
.476705575+ 04 
.7* 9 06 8 4 35- 4 
.45 n 588 075+ *4 
.74 ‘'"166735 - r 4 
.46 2559 035+ 04 
.74*840575-' 4 
.4757*5 ^85+ 9* 
•18121° 665- 03 

• *4 7 60532E + 0 4 
.7 4" 368 5 05-0 4 
.461=58 135+ C4 
.74 5 1205=E- 3 4 


T 

.166 "'497 ' 5 + " q 
.524 076925 + u4 
.1673622*5+ 09 
.519 78 79 2 5-+ ’’* 
.166 1962*5+ -9 
•550359465+04 
.166 '81865+09 
.517866775+ V4 
.1674356 :£+r9 
•5172=4475+04 
•1661 78265+09 

.5491:9695+04 

.166 983945+09 
.511607255 + 04 
•167513725+:9 
.514388385+ 04 

• 16620;1655 + ,! 9 
.546884335+24 
.166 u553 35+ u9 

• 5 u5 2 8 6185 + 14 
.1675836 "5+ 09 
.512714225+:* 
.1662589 15+ v9 
.544 *93995+ n 4 

• 165=97525 + <-9 
*499896265+ 04 
.167647385+69 
.510731625+04 
.166317715+9 
.542114135+' 4 
.165=14625+09 
.492435455+ J4 
.167705615+09 

• 5 03 £9 3 7 2E + 2 4 
.166 328 1 45+ '9 
.539864415 + f; 4 
.165*129 1E+99 
.485° "1395+ 04 

• 167761415 + "9 
. 5 "715 0325 + : 4 
. 166251 O0E + ' 9 
.537*437 .5+04 
.165711 245+ 29 
.479 2 029 35 + 04 
•1£7R15265+ P 9 

• 5 0 5*376 '.5+14 
. 166*06445 + C q 
.53485*275+^4 
•165591785+-9 

• 472619 5 55 + '‘4 
.5 67*718=5 + ':° 

• 5 " 3 7 1 3 3 9 5 ♦ C * 


M 

. 125395375+05 
. 263514965+01 
• 121 652675+25 
. 267377 = 15+11 
. 1299845*5 + C 5 
. 248281785+01 
. 115982535+05 
.2674649 0 E+O 1 
. 957478665+04 
. 269452165+01 
. 124412565+05 
. 249846725+31 
. 1 172 2 ? 345 + C 5 
. 271419355 +C 1 
• 912824445+04 
• 27143245 E +01 
. 119022165+05 
. 251642335+01 
. 990245645+04 
.2752878 35+01 
.* 72226225 + 0 * 
• 273277915+11 
. 113602425+05 
. 2535 = 03 1 E+C 1 
. 913928845+04 
. 279382955+01 
.834 = 414 8 5 + 0 4 
• 275 C - 13865 + 01 
. 108293755+15 
. 255540935 +J 1 
. 6428 * 1435+04 
• 283411335 +ul 
. 802328535+04 
. 276669525+01 
. 10335 = 225+05 
. 2573494 15+01 
• 776688985+04 
. 2974 = 4145+01 
. 767661 : 65+04 
. 278278495+01 
. 984783415+04 
. 25927 = 265+01 
. 715263095+04 
. 2916 * 425 E+C 1 
. 736343355+04 
. 279978255+01 
. 935241375+04 
. 2613628 . 25+01 
. 653172 * 65 + 0 * 
. 295888315+01 
. 7 057592 75 + 14 
. 28 1 * 99695*01 




A- 2 2 


ITERATION NUMBER 2000’; DTIME = . 24 76 7 2 - 2 E-S & TI W E = . 346 7359 IE-92 


NODE 

X 

Y 

R H 3 

U 

V 

SOS 

r 

P 

13 



GAM 

cv 

T u X Y 

THX2 

M-J 

T 

U 


2761 

.1576337 3E+ 01 

•1124779 2 E* ''l 

• I5363237E-r«l 

•11352523E 4 25 

• 5 23878 77 E 4 04 

.474894 7 6 E 4 84 

. 16642«5*E+T9 

.eS9267o2E+04 

*? 



•12119891E+01 

. 6 6 1 9 6 5 a 7 E 4 0 G 

.247716528+ C 2 


• 1 8 6 6 3 7 5 2 E - 0 3 

.532 5 39128 + -J4 

.2632791 18*01 


27 73 

.169 09696E+01 

• • 

. 1 l a 53265E- :l 

• 133*7 25 3E 4 05 

. r. 

• 44458751E 4 .4 

.16-5*666 i E 4 v 9 

.6 152595 " E + c* 4 

1 



. 12 1326 14 E 4 r -l 

• 6548555 ‘ E 4 C 3 

• ■ 

. 9 

.74 *‘'172228“ *4 

.46587956E 4 "* 

.30 0216568 + 1 1 


2765 

•16529221E+G1 

.74 j 1 6 7 46 E* ,0 

.1238 193 C E-31 

•12338 C5oE 4 75 

.4 2935599 E 4 24 

• 46 1313 39E 4 84 

.16792595E 4 09 

.675464138+04 

9 



.12125032E+01 

.65 72 12C3E+ 0 0 

•19187559E 4 C2 

. 

• 7454 04 Q 4 E- 2 4 

.5ri90076E + t)4 

.283183428+01 


2797 

.1631S8U0E+C1 

. 11495155E* U1 

. 14694451E-01 

. 1143 £ ^ 8 9 E 4 5 

♦ 5 2 495021 E 4 04 

.47413543E+ 04 

. 16641295E 4 09 

. 64 7 6722 - E 4 C4 

9 



• 12 12 OB 39 E 4 01 

• 6565 2254E+- 2 9 

•244£567?r*0? 

* 

. 

.19239979E-'3 

.53 >35169E+(.4 

• 265C463 : -E + 01 


2 8 09 

• 1745645 IE 4 G1 

. c. 

• 111425198-01 

• 1 345 1 715E 4 C5 

. ? 

.44141975E+ 34 

.16533163E+''-9 

. 555696538*04 

1 



• 12 1435 uSE 4 01 

• 65149828E 4 ■. 3 


. r 

.740070795-24 

.458671558+ 24 

• 3 2473751E 4, jl 


23 21 

.17 394955E+ Jl 

•75663647E 4 ^ 8 

•11873949E- 1 

• 12432634" 4 75 

.42 493913E 4 04 

.4 6 26 795 CE 4 14 

.16799299E+ 39 

•64517 2 928*04 

9 



• 121363 03 E 4 " 1 

. 65371474^ + '•) 

. 187T f 731E 4 ?2 

. c 

.745597 37E- 24 

.500 05 C 15E 4 14 

• 2849 234 7E 4 C 1 


2833 

•1&694 8758+01 

• 117517 C 3 E -4 01 

.14 02984 IE- Q1 

. 1152469?E 4i :5 

.51919752E 4 14 

•473343 49E+ C4 

.1664430 5 E+ 09 

• 8 052981 5 E+J4 

9 



. 12132219E+ 31 

• 65494698E 4 0 0 

. 2421 o 219E 4 C2 

n 

.icg. nqg 26 E- 03 

.5241 O59SE+04 

.266953698+01 


23 45 

•13Q31456E+C1 

• J 

.1037752 BE- 31 

• 1 355 7 1528 4 C5 

n 

•43785458E+C4 

•16520266E+09 

.5088561 vE*04 

1 



•12152166E 4 01 

•64886279E 4 30 

• ‘J 

.0 

.740067858-^4 

.4511779 :E+04 

.3 096268 2E + U1 


2357 

• 176 9 9 71 4 E 4 01 

.773861728+ 33 

•11356269E-01 

• 1 25 2 c 32 IE 4 05 

.41 27 a 573E+ 34 

• 4 598 1 7 23 E 4 04 

.163 0378 5 E + 09 

•6 14334968*04 




. 1214775 JE 4 01 

• 6 5 0 2 0 1 6 3 E 4 ) 0 

• 13236538E + 92 

. 4 

.745517368-44 

•49775211E 4 -4 

• 286 = 718 SE + O 1 


29 69 

• 175 0 071 4 E 4 Cl 

.12 ;1 3 4 43 E 4 31 

. 13223529E- •;! 

•116167'1E 4 05 

.51 73654SE 4 34 

.4 72250298+ ^4 

.1664949 r Z+ -9 

.76 L4715*E 4 Q4 

9 



•1214439 3 E+Ol 

• 65122799E 4 - 0 

.24 0 96447E 4 f2 

.•J 

• 1° 3468 428- 03 

.525 17263E+ ;4 

.2692788 OE + O 1 


2381 

•13633 09 3E+ 01 

• 

.96*S8?17E-' , 2 

• 12£6 7829 r4! 05 


• 43 397371E 4 04 

•165 2533 1E+59 

•46447548E+04 

1 



. 1216uG 04 E 4 .1 

.64c,4=>Q65E 4 . ? 

• J 

. r 

• 74 C 063 4 OE- 04 

.442°3552E+ '4 

. 3 146538 4E 4 21 


29 93 

.13317297E+:i 

. 79176579E 4 O': 

. 1 J 8 3 7 * Q 6 E - ~ 1 

. 12625 Q 73E+C5 

.4 r 5 1 5 947 E 4 04 

.4 585 01 8 6 E 4 04 

.168 13412E+ r 9 

•58253232 E* 04 

3 



•12155£59E 4 31 

.647747168+00 

.177 a 1 1 06E +02 

. ! 

.745291948-04 

• 49458 499E! 4 -'4 

.28920524E+01 


29 0 5 

.1913989 1E+01 

.122965 G1E + 31 

.126322C5E-01 

. 1 172784 1E+05 

.51645131 r+ 0 4 

.4 7 0792 78E+ 04 

• 16655612E 4 '2 9 

•7160278 c E+04 




. 12153567E 4 u 1 

.648439168+ r- 

. 238 ’3 0328+ C2 


• 19 7834 ? OE- 0 3 

.52155313E+04 

.2715036 " E + O 1 


2917 

•192952S9E+01 

. 0 

.89445949E- '2 

. 1376921 08 + 05 

. " 

.429921 31E 4 34 

.16*6414 5 E 4 -9 

•4222762 8 E 4 0 4 

1 



.121694 95 E+ 31 

• 64395891E+ O' 

• • 

• 2 

.741185658-04 

• 434*'0 48?E 4 J4 

.3202728 . E 4 0 1 


29 29 

•13977 23 3E+G1 

.839646938+ 06 

.1 332119 8 8-01 

•12722*768+05 

• 3 C 798517 E 4 94 

.457162 q 3E 4rt 4 

.16922753E+:9 

.551159a*E + 0 * 

9 



. 12164383E 4 01 

• 64516677E 4 -. 3 

•17366 65 2E 4 C2 

. 2 

.755918998 - *4 

.491 36299E+(: 4 

•291594 8 5E* 01 


2941 

•19831032E+01 

.12577341E+21 

. 12 J58714E- ■: 1 

. 117 c 2538E + f’5 

• 5121 ;759 E 4 24 

.469661 348+ 34 

• 16£54113E 4 '-9 

.6797608 OE + L'R 




•121621398+01 

. 6 4 5 6 6 7 B 4 E 4 0 0 

.234719118*52 

m r 

•416909948 - 5 3 

•5196914 6E+24 

.2737593 2E*u 1 


29 5 3 

•199525278+ul 

* 

• 8 3 1 3 9 9 6 9 E - ' 2 

. 1 3 86 7 5 75E 4 25 

# 0 

.4259151 78+ 04 

• 16458541E 4 '-'9 

. 394556058+ 04 

1 



•121769S8E+G1 

• 64143769E + 3 0 

. C 


.740287 0&E- ^4 

.426 0554 1 E+ u4 

.325594768+01 


2965 

.196369508+01 

.8277 „924£+ i- 

.9&36469 'E--.2 

. 12 614 369E 4 05 

.3912253 3E 4 34 

.455862998+ r 4 

.16332255E+T? 

• 521925*5E 4, ;4 

' 



• 121729 09 E 4 A 

• 6 4 2 6 4 6 1 6 E 4 ■' 0 

• 16977537E + 32 

. 5 

.765239378-84 

• 48823744E 4r: 4 

. 29391 01*E 4 01 


29 77 

.1?‘+62265E 4 91 

• 12= : 5472 r| E 4 01 

• 1 1 6 2 775 9 E- "'l 

. 1 1 87 3 3 14 E 4 05 

.5 r ’ 3 7 3 5 3 9 E 4 34 

•*69063678+ 34 

• 16646583E 4 *- ; 9 

• 65 33421 48 + ^4 

9 



. 1 2 1 7 ' 6 5 3 E 4 '1 

• 64331 6 45E 4 v ' 

. 22989471"+ *2 

. 1 

.62* '91 18E- ’’ 3 

. 5 1 7 1 7 0 3 E 4 ~ 4 

.2749794 'E+. l 


29 89 

.2 :6 G9947E+ A 

. 0 

•77*926 C 58- 12 

•13959 Q 03E+05 

. r 

.421955658 + =•* 

.16429594E+-.:9 

.351921738+04 

1 



•12135483E 4 31 

.63 8 a 3544E 4 ’: r '- 

© 

. 2 

.74 0277928- '4 

.417885^18+ ’4 

•33083617 E+Cl 


30 31 

• 2.29683 0 E 4 & 1 

• 8 4525213E 4, 3;j 

.9333491 0E-U2 

.12942 ' 9 6 E 4 0 5 

.3 8 5 1 2 5 3 7 E 4 n 4 

.454613 '48+ '’4 

. 1 6 0 4 1 3 8 4 E 4 n 9 

.4 Q 4817J*E+"4 

9 



. 121 8 14 3 7E 4 _*1 

.64 1 2479E 4 99 

. 1£61539*E 4 02 

• Z 

.765 01 0238- '4 

.48523 0 6 6 E 4 " 4 

.29617:;04E + 'l 


3313 

.2 012366 28+ 31 

•13128724E 4 01 

. 11293854E-'’i 

• 11 q 5"944E + /1 5 

.4C284293E 4 24 

.46°7 c, a 26E 4 ^4 

.16925937E+'? 

.6333963 *E 4 C 4 

9 



•12179199E+ 01 

.64 178*458+ 0 0 

. 2241 067SE 4 2? 

* j 

•628987528-03 

.516 1 5 9 5 1 E 4 1 4 

. 2757582 1E+G1 


3125 

.2126761 0E+G1 

. 9 

.7243425 E - ' 2 

. 1 4 04 7 24 8 E 4 ' 5 

. 

.4 1 8 047 948+ r 4 

»163 n 8655E 4 ‘ ? 

.322659268*04 

1 



• 12 19 35 93E+ u 1 

• 63£4516££ 4, ‘‘3 

• r 

. 7 

.740263948-04 

• 4 : . a 89944E + ‘ ; 4 

. 336 ::2'.59E + o 1 


3037 

*2 J957C55E+G1 

.86 273 3 07 E 4 0 !. 

.8962431 OE- 02 

. 1299586 or 4 05 

.379 25593E 4 ') 4 

.4 5 3424 T 08+ ^ 4 

.16351*398+29 

•46931599E 4 ;4 




. 12 199968E+ 01 

•S3762213E 4 "0 

. 16 27 24718+02 


•76*270 a 9E-'4 

.482 365*6E + '.4 

.2983471 E+Gl 


3 0 49 
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3 .i0145l7E+"4 

.4165 06865 + 01 


1 75 <05235+ .-9 

. 129961655* l- 4 

9 

4126661 7E + 1 4 

.3554268 15+01 


170 1414 :e+ : Q 

• 114711.i 55 + 9 4 

9 

3 8 C 3 C 7 8 9 E + 0 4 

*37 C35 G 465 + 91 


1712701 CE+09 

.132222.75+04 

a 

392 57768E+ 04 

• 3633 '-5585 + 21 


1 7 [ 628295+09 

. 13 29 20 8 3 5 + 04 

9 

394 19794E+: 4 

•3611 Jj9 55+91 


170 46S32E+C9 

. 159175695+C4 

9 

4 56 54278E + .. 4 

.352261285+01 


1 7 0 1 5 3 2 7 E ♦ 9 

. 17732S8 35 + 04 

9 

418 5623 * E ♦ 0 4 

• 34 36979 <5+91 


17" r 39515+ '9 

.2'73£2325+04 

9 

427782195+ * 

• 3275221 c E + , j1 


170 .'26465+ ’ 9 

.2124451 55+ 0 4 

9 

432 '0967E+ 4 

•234802655+51 


16987*3 E + r.9 

. 22*3562*5+04 

9 

435174235+ '4 

• 3325645 2 5 + 1 1 


169811065+09 

.23: 391-C6E + 04 

0 

439 2 9 9 2 5+1 4 

.33:943935+01 




A-25 


ITERATION NUMBER 2 CO O': "TIME = • 24 76 72 ' 2 E- 3 & T I M E - • 2 *6 73 5 q 1 5- " 2 


NODE 

X 

r 

RH3 

'J 



GAM 

cv 

thxy 

3 5 76 

.3 j 9 6 1 7 * 1 E + 0 1 

. 10 37 6 5945 ♦ ,1 

.5 -.65 7512"- .2 

.139=220=5 +05 



. 1228 7 755 5+ ; 1 

• 6 1 ■- 3 J 3 9 4 £ + ■ 

. ll c 7769°E+f-2 

35 77 

•3j9*5=:iE+Gl 

.11 ‘'19 , 55 E ♦ 0 1 

.511 !7: 5 IE- '2 

• 1391 68445 + ^5 



. 122576 23 E+ 01 

.61 3 3 3 9 1 1 E + -0 

. 12539 fl 1 5 E + 0 ? 

35 76 

.3 9307735+Jl 

. 11632331E+ 01 

.525 = 0413E- 32 

. 138* •* ! S623 + C5 



•12 2874=75+ -1 

.61 03727 :e+ ' : 

.132149885+02 

3579 

• 3- 91628 7E + 21 

. 12219 6 02 E+ Jl 

•537*9=385- 02 

. 137647425 + 135 



. 122873 775 + 0 1 

• 61 34 C 4 6 6 E + 0 0 

. 138473345+C2 

35 8 ; 

• 319 0233 05 + 21 

.12733*235+ Jl 

• 555715115 - "2 

•136 895765+0 5 



•122872615+01 

.61 34 357 7 E+ 0 0 

. 142652675+C2 

3581 

.3 J986994E + Cl 

• 13 32 & 363E ♦ 01 

.57*266325-02 

.1361458=5+05 



•1223715 3 E+ul 

.61 04655 "E + 0 5 

.146335455+02 

3582 

• 3 08 76 08 1 5 ♦ Cl 

• 133*9.5 72 5 + 01 

.55333U45-C 2 

.13544919B+C5 



• 12287 .) 4 3 E ♦ 0 1 

.61 '*9*195+ Of- 

. 149368l7 r + (32 

3583 

.3,9635955+01 

. 14355 782E + Jl 

.695747*25- 02 

. ! 2 47 7 626 E + ?5 



• 1228 6 9 39 E + 0 1 

.61 -521935+ 00 

•154 P 2 7925+92 

3584 

• 3 0 3 5 1 4 9 4 E + 01 

.1*8*635*5+ 21 

• 608 3 89*55 - '-2 

. 13*107395+05 



•12 2868 33 E+01 

.61 05 488 1 E ♦ 3 D 

.16151*435+02 

3585 

.3 =8397425 + 01 

.153227=55+01 

.612287365-02 

. 133362705+05 



• 12286 7 41 E* 21 

• 61 057493E+ m!0 

. 168743 345+C2 

3586 

•3932830 45+01 

• 157B64B.JE+ Jl 

• 61567672E- '-2 

.132521 :55 + 85 



.122366*65+ 01 

.61 06C035E+ ;, r: 

. 17365 7 865 +?2 

3587 

.3 J&1715GE+C1 

.162386705+ 01 

•62975213E-02 

.131497585+05 



• 122865535+ 'J 1 

• 61 0 625145+ 0 ' 

•1776369^5+02 

35 88 

*303062515+ Cl 

•1668 05 255+91 

.636 = 48675- "2 

.133220185+05 



• 122S6462E + J 1 

.61 ‘‘'6 49365+ OC 

•18 ?5 18925+02 

358? 

.3 .-79556 35 + 01 

.17113 1 19E+ 01 

.65 09 26= 2 E- *2 

.128503775+05 



• 1228637*5+ 31 

• 61 r, 673J85+0 0 

.184752115+02 

3 5 9*. 

.3 J735114E + G1 

.17537*515+ 31 

.65 32333 IE- "2 

•12617 05 4E+ 05 



• 12286237E+ = 1 

.61 0696355+ OC 

. 188301565+02 

35 91 

.3 J 77*8255+ 01 

• 17°54454E+ 01 

•662123 8 2E- 3 2 

• 12 2 8? =5? E +05 



.122862 J1E+ 01 

.6 1 ’7 1 9 2 25 ♦ ? 

. 193165385*02 

3592 

.3 7547015 ♦Cl 

• 18365 0 05 5 + 3 1 

.655727645- r 2 

• 118 327115 + 95 



• 12286 1 17E+ C 1 

.611741735+ 3 

. 19 6*265 0 5 ♦•■'■2 

35 93 

• 3 :75 4 71 2 £>01 

• lc 76 9 9 32 E+ !J 1 

.66=3063*5-2 

•1119?1125+ r 5 



.12286 03* E + ol 

•61 0763=*E+ T 

• 2 ' ; 1 4 = 1 8 5 r + 02 

35 94 

.3 . 74 464 35+01 

.19170 '235+01 

♦659333295-02 

.103127445+05 



. 12255952E+ 01 

.61 37859 = 5+ " 

.204663785+02 

35 9 5 

.3 = 7 3 5 *;7 5 5 ♦ 01 

.19566 ?29E+ 01 

.687586575- 32 

.81 1676215 + 04 



. 12285&71E+ 0 1 

.6 1 9 0 7 6 1 E + -‘ 3 

.21 0C1 4715 + 92 

35 96 

.3 07 2 5 39 05 + 61 

.199586 74 E+ 01 

.693 '.4595E- 02 

. 76 = 3 0 7075 + 04 



•122857905+01 

.61 ''3 291 55 + •: 

.2132783*5+02 

35 9 7 

• 3 0 71 5 771 5 + 01 

•203*86585+ U 

* 74 24 39 4 7 E-*. 2 

.6360397*5+0* 



. 122357105+01 

.61 .‘850555+ C 3 

. 2198 = 3*35+02 

3599 

• 3 J7 362005+81 

•207366605+01 

.74 33584 E — 12 

• 5153 1 ^ = 65 + ** 



•122356315+01 

.61 :8 718-4 5+ ‘ ? 

.22 07 2«3 8 5+02 

3599 

• 3 0 6 = 66615 + 31 

•2112335 3 E+ 01 

.56 J 29162E- 32 

.35287:345+0* 



• 122855525+ ->1 

. 61 0 39335E+ 0 ' 

. 22378 6475+ 0? 

35 :: 

.306371595+01 

.215 0-93*65+1 

.1.61393 * E - - 1 

• - 



• 12265472E + 01 

.61914245 + - 


36 01 

•3169:5295+01 

• 0 

• 3 2 5 ' 8 4 9 4 E - 2 

• 1 4 °9 5 £94 r + "5 



. 12295473E+ 0 1 

•6.82*8725+ 0 

• • 


V 

sns 

E 

P 

IB 

TMXZ 

MU 

T 

w 


286681255+04 

.436231 = 05+'’* 

. 1696 = 38 3- E+‘ = 

.2*3837855+0* 

9 


•753 =03135-0* 

,**2=726 8 5+ T * 

•327*12355+01 


3 "9523375 + ff* 

.* 3 7288 735+0* 

.169656965+6= 

•2*71977 1E+0* 

9 

n 

.7587*2175- 0* 

.**5 126335+ J* 

.326 =29295 + 01 


325 -'16*9 E+ 04 

.*33720 07E+ 0* 

.1695**515+== 

.2560*31 1 E + 0 * 

9 

r-. 

.767679735-"* 

.**8 :*96*5+:* 

.32*06*275+01 


3 3-93 10505+ 0* 

.** =276 = 25+ 0* 

.169*59*75+:? 

.263551555+0* 

9 


.78 =869 8 IE- 0* 

.*51 239595+r* 

.321=96515+01 


3*8059275+0* 

.**2123275+ :* 

.169 31*095+ 

.27*782225+0* 

a 

r 

.8*1516915-?* 

.*55:*67:e+c* 

.319*80155+01 


355*92605+0* 

•***16573E+ C* 

•169202515+0? 

• 28653 ?6 =5 + (J* 

a 


.833**8315-0* 

•*5=25**85+ "* 

•3168-38=5+01 


3613345SE+ "4 

.**6 058 = 25+'=* 

.16= =65375 + 09 

.298652185+0* 

9 

r* 

.8835**28E-'i* 

• *631818*5 + '•* 

.31*276835+01 


3 7335709 E* 04 

.4*727895^. ,14 

.169=69*75+ = 

•3 265*88 8 5+0* 

9 

7 

.960862595-0* 

.*657229*5+0* 

.312672935+01 


3 9 8 39 584 E ♦ ^4 

.**7879685+ 0* 

.169 8577*5+09 

. 3C8716G55+G* 

9 

n 

.1 079713 0E- = 3 

.*66 978585 + * 

.311731385+01 


4 C453443E+ ^4 

• **8*85=25+=* 

.16969=085+09 

.311*97*15+0* 

9 

r 

.127357 =25-03 

.*682*73 ^5+ ( * 

.3107*1*55+01 


41442697^04 

.* *9516225+0* 

.1633190=5+ j9 

•31*7 J 69*5+0* 

9 

n 

•15=6185 05- p 3 

.*70*9* 795 + ''-* 

.3 j 83876*5+01 


4214 9004E+04 

.*51073675+0* 

.167 76 11 65 +09 

.32413766E+C* 

9 

* 

.21 7*8885E-"3 

.*73672625+''-* 

.30612*555+01 


4 2 4 4 1 529 E* 04 

.*52758865+ ** 

•166762635+ “9 

. 330296385+0* 

9 

r , 

• 31 = 3 = 0 035- •. 3 

.*77222975+ C* 

•3 025 05325+01 


42934945 r + ->4 

.*5*5*5 225+ *4 

.1653332*5+0? 

.3* 021S28E + 0* 

9 


•51*05*715- 03 

.*60999615+0* 

•29837 C75E+C1 


4 3 1 2 6 0 7 ? E ♦ 04 

.*56*31*65+5* 

.163 163 85 + 0 9 

.3*268*665+0* 

9 

?: 

.88118686E-=3 

.*85 :033*5+0* 

.2=2059395+01 


43 075144 r + 04 

•*567?5=05+C* 

.1596=6 035+0 9 

• 352**30 0 5+?* 

9 

r\ 

•166966635-62 

.*89852* ■* 5 ♦ C * 

.283886175+01 


42233665E+ *4 

.*6 1 276 *25+ =* 

• 15* 6796 2 5+ :9 

.35296366 5+ u* 

9 

*" 

.3 3558* 59E-02 

.*95361275+0* 

• 272371 r ‘ 15 + 01 


4! 1662° 1E+ j 4 

.*63902*55+=* 

.1*76868 55+ -.9 

. 36*3856 *E + 3* 

9 


•7638*6 1*5-02 

.5 01 0238 25 + o* 

• 2569379 £: E*C1 


384899 36E+ 04 

.*65622*75+ 0* 

.137778775+09 

.3 61625* 5 E + 0* 

9 

n 

•186198235-01 

.50 *7*62*5 + C * 

. 236* =5** E + 01 


3 4 °? 9 657 E+ J4 

.*6*521 *35+ =-* 

.12*516295+0? 

. 3753*28 jE + Q* 

9 

1 

.5286*2 !25- r 'l 

•50236536E+:* 

.210226*75+01 


i 73 71 69 E* 04 

.*57722 1 1E+ =* 

. 1 08 7 0 72 . 5 + v 9 

•367331? 75+0* 

9 


• 1*29 7=535 + r -0 

•*8776 n 66E+0* 

• 18? *29815 + 0 1 


25674875 E + 04 

.**6679 175+ 5* 

.945 7*28 ?5* 08 

.37*75*735+0* 

9 

7 

• 2 = 663=3 65+0-0 

.*6*52=835 + '* 

.153556665+01 


~896.31SE* 04 

• * 3 **38 65E+ '* 

• 8267376 = 5+ '' 8 

.355176865+0* 

a 

7 

.33=919965+ 0 

.*39*126 5* * 

• 127997625 + -J1 


4 9 ? 0353 E+ 14 

.*1*122915+=* 

. 68 * 1 0791E+ ' = 

.37672*7=5+0* 

a 


.21 = * = 6 3*5+ 00 

.3=926:; 115+ * 

.92*8*89 = 5 + C-: 



.35 =55 0 125+ : * 

.*6 0 *16675+ : 9 

»3*713*S2E+u* 

2 


.7* "TCO -?r- ’* 

.3 ::=3235E+ * 

* 



.35516? 125+ :* 

.157127535+. ■? 

.1:3656115+0* 

1 


.7* 0591 * “5- :■* 

.2 = 3**:-3F5+r* 

.*22223*95+01 




A- 2 6 


ITERATION \l JM358 2 0001 7! T I M 5 - . 2 4 7 6 7 2 ■' 2 E - 0 6 TIME = . 3*S 7 3 59 1 5-52 


NODE 

X 

Y 

RHO 

U 

V 

SOS 

E 

P 

15 



GAM 

cv 

THXY 

THX2 

MU 

T 

M 


3 6 1 3 

.2 15 2 134 6 5 + >11 

.11 1652616+ U 

. 4 9 4 = 6 v C 4 r - . , 2 

♦13964731E+C5 

• 3 r 644914E + 04 

.436318535+04 

.16972612E+C9 

.238 0 42795 + 04 

9 



•122932365+0 1 

• 6 G c 9 4 2 9 * E ♦ - 

•12 3 77*325+02 


.8068617*5-04 

• 442 a 52565+ 04 

• 32766 49 95 + 6 i 


3625 

.314724755+01 

. 173* "339 5 + U 1 

• 63279618E--'2 

. 1269 ‘7 7*55 + 25 

•425674315+04 

.453317905+ P4 

.165 2 4 76 f .E + 9 q 

.329533465+04 

a 



•122919396+71 

•679172L7E+' 9 

. 18 2T4 34 8E +f 2 

# 'j 

.16*672655- *2 

.4792*3325+04 

.259136195+ 01 


3627 

• 32565 03 05 + 01 

. 0 

•31278952E-22 

• 15C2S3956+"5 

.0 

.351153405+04 

•1564894GE+ 9 

• “745355 36+G3 

i 



.123 01 9 7 9 E + , 1 

• 6C6764 8 2E + . 0 

. n 

r 

.74 7572 015- 74 

.2«672717E+04 

.4279723 75 + 01 


36 4 9 

.322984746+01 

.11217 2 2 06+ Ol 

.479:88135- 2 

• 14 91025 OE +05 

.3 -I 3242325 + 04 

.435404855+ 04 

.1697928 2E + 09 

•229525996+C4 

9 



.122988 645+ vl 

.6 7.73 5 02 86 + -7 7 

• 12212892E +02 

"i 

.8 05112 445- 04 

• 44 " 899 1 7E+ 94 

•3292261 95+61 


3661 

•32151C56E+U1 

.175653135+ 21 

.615514575- 92 

. 12°2 C 6315 +C5 

.422 006955+ ”*4 

.453118 *25+04 

.16515628E+ '9 

.319477065+04 

9 



• 12297 633E + 7 1 

.6:276 74 7:5-*- "0 

•18076 *355+52 

r 

.161187Q75- 72 

.47755u24E+74 

. 3r , 016203E + 01 


3673 

• 3324 129 7E+ 01 

. G 

• 301313925-'. 1 2 

•155570995+95 

r. 

.347152235+04 

.155811025+09 

.917791045+03 

1 



.123066 99 E+0 1 

.6:5284295+99 

• C 

.0 

*7 4 0557 635-0 4 

.280173265+C4 

.4337319 ,.5 + ei 


36 85 

•329773756+ 01 

. 114536626+ 01 

.4647 7516E- 72 

. 14 054*225 + 65 

•2 C 9799295+ 04 

.4 34559535+ 74 

• 1 69 8546 25 + :, 9 

.2216699 45 + 74 

9 



•123945066+71 

• 6 058 6122E + n 

. 12 0*14265 +09 

.0 

*803917 635-74 

.438988975+ '4 

.337693895+71 


3697 

.326314145+31 

• 17788 085E + : 1 

.59 °R 53155- j 2 

. 12°669165 +C5 

.418560205+04 

• 4 5 2 429 57 5+ 04 

• 165 l'6239E + .’i9 

.309667685+04 

9 



• 123 032935+ 0 1 

• £ t< 6 1 8 f 7 7 6 + : 9 

. 178896 75E + 02 

• 0 

.158823 775- 02 

•475882255+74 

.301 1677 1E+01 


37 "9 

•339194216+31 

• J 

.29 989464F- : 2 

.150845435+05 

# " 

.347236105+04 

.155152135+09 

•86512 017E+C3 

1 



.122123356+01 

.6038 069 15+ " 9 

• . 

. r 

.747543685- *4 

•2736957 15+0 4 

.43948^1 GE+C1 


37 21 

.3365814 06 + 31 

•115957156+ -1 

.45141273E-02 

. 14 3967335 + 05 

.296043415+04 

• 4337 Q 0 3°E+ r '-4 

.1699297 7 5+ 09 

. 2 1 4 46 q 3 76 + 0 4 

9 



.1231 01646+C1 

• 6 n 4 3 752 9 6 ♦ ' 9 

• 1 1 8 6 0 249E + 02 

.0 

.8 02706 025- 

.437236? 15+ 04 

. 332 -7554 1E+C 1 


3733 

• 3 3 5 1363 9 6 + 01 

.18 .9 0S 7 JOE* 0 1 

• 583 1 530° E - 2 

• 1 3 0 C 2 58 36 + °5 

.4 14842865+04 

•451 8 33745+ r 4 

*164966935+99 

• 3 O 0 5 7 6 9 7 6+ 0 4 

9 



•123789636+ 21 

• 60469fiC' Q E+'.‘9 

.176 Q 51.UE + 02 

. r 

.156317165-72 

•474349725+ J 4 

. 3 J2 08 51*5 + 01 


3745 

.345994936+01 

. 3 

.26166961E-";2 

• 1 5 1 1 3 47 q E + 95 

.0 

•33 G 474815+04 

•15456986E+79 

.819 04976E + 03 

1 



• 12 2 1 7 9 87E ♦ C 1 

.6 C23326 1E+ 7 

. 1 

. " 

.74 0532725-"4 

.267 607525 + C4 

.*45271796+01 


3757 

• 3 434 686 7E ♦ 31 

• 11 7364 .,5 5 + n 

.439173675-72 

•14137167E+05 

•291938555+ 04 

• 4 3 3 ' Q 7 4 05+ 04 

• 169 Q 56B :'5+i/9 

.207892325+04 

9 



•123158336+ 91 

• 622892355+' 9 

• 11667815° + 2 

# r* 

• 8 0 1734 115-04 

.435641625+74 

.3333773 15+01 


3765 

.341978246+01 

.18227159E+71 

.5685661 IE- 92 

.130367125+05 

.4 1 '*91456 5+ '* 4 

.451214 175+ 74 

.164 °46775+79 

.2921539 16+C4 

G 



.123146496+01 

.673292465+ .9 

. 17494747E+G2 


.154328515-02 

.472895595+ "4 

• 3 0293775E + 0 1 


37 6 1 

.3529 16356+01 

• 3 

.273587945-72 

.l£145466E + ' : 5 

. 7 

.335904045+ 04 

.154 094655+79 

.77854333E+J3 

1 



•123236566+01 

.6276 6 0885 + u " 

• - 


.74 .7523 785- *4 

.261887885+04 

.4578867 :E+G1 


3793 

.35 02563 uE+Ol 

• 1 1 S75 7 6 36 + 91 

.427897555- 02 

. 1 4 1 75 8255 ♦ 05 

.28 7492195+ 14 

.4 32478 2 05+ 74 

•169996365+09 

.201882*65+04 

9 



•1232152SE+91 

•601412295+^9 

.114643585 

. .0 

•3 G 0802765-74 

.434 197 3 5E+ 4 

.3344541 26+01 


33 C5 

•343B4LS3E+01 

• 18*436235+ 9 1 

.55471C75E- "2 

• 13 369362E ♦ 75 

.4 ‘6718045+ 04 

.45 07 71 295+ *4 

.164731075+79 

.264267176+04 

2 



•1232 9 3 52 E + 01 

.671717695+"? 

. 172861225*02 

. 0 

. 15 23 02 155” n 2 

. 471624665+ n 4 

•333695436+01 


3317 

.359658536+31 

. 3 

.2664 08325-02 

• 1518 t <07 Jf + 05 

, 2 

.332494855+74 

•15371157E+C9 

.742459346+03 

1 



• 12 329 2 42 E + 01 

•399391855+ r 9 

. r 

, r 

.749517635-74 

.25648148E+-:4 

.45655.7516 + 01 


33 25 

.357125425+01 

.12 213 8 94 E + 71 

.417482655 - ‘*2 

. 14212878 E + 7-5 

.282741615+ 04 

.431922 035+ *4 

,170728686+0° 

•1 96371456+ c* 

9 



• 123272 3 7E+ ’J 1 

.599934655+ 9 

.11251 15 Q E+:2 


.8 0 0283795” 0 4 

.432882925*94 

.335579246+01 


3341 

• 3557 244 95 + (Jl 

.16658 ; 12 E + 91 

.542277525- 2 

.131 007395+05 

• 4:2 016485+04 

.45 "234825+ ■’*4 

•16459728E+C9 

.277185166+04 

9 



• 12 3 26 7 73 E+ 7 1 

.6 .<923 55 85+ 9 

. 17 0 594*45 + 7 2 

n 

.151356845- 72 

•47 G 412236+04 

. 3 7-4 36 75 6E+ 01 


3953 

.3655233 05 + Jl 

• J 

.25577733 5- ’2 

.152148535+05 


• 3 2 q 1 °5 3 35+ '4 

.153 37 7465+ *9 

.779353466+03 

X 



.12335 7485 + Cl 

.597925225+ ‘0 

Q 

, 2 

.74 "512 :<4E- 74 

.2513 0C955 + 0 4 

.*621833 ! E + 7 1 


3365 

•364716936+01 

•12150562E+71 

.4 078 22 = 1 E- "2 

.142 43 440E-*n5 

.2777516 05+ ^4 

.431426485+04 

• 17 2 ;i56445+ :9 

. 191298985+04 

9 



•12332 ° 65 E+Ci 

.592459P6E + _• 'J 

.11 73^6075+02 

. 0 

.79°737 115-04 

•4316 c 11 c 6+04 

.336479886+01 


39 77 

.35263 77 86 + 31 

.1887 94 04 5 + 01 

.52635642E- ?2 

•13131 23 55 +05 

.39757462 5+ ' | 4 

.44 ° 7°9 965 + 74 

.164 476725+75 

.277187385+04 

3 



• 123316 13E+' ! 1 

.59=755965+ 0 

. 1694 47475 + "2 

r. 

.15 7 076 76 5- 0 2 

• 4692 8 7165 + 0 4 

.375 7,22445 + 0 1 


35 69 

•373411326+01 

. 0 

. 25335°945- .2 

. 152467465+75 

• 0 

.3 2588 7585+ "4 

• 1 52 °9 5 0 1 5 ♦ t'9 

•67765 05 6E+C3 

1 



. 1234 77725+ 0 1 

.5 964 6 7 915+ . 7 

• 1 


.74 75 09 3 05- r 4 

.246 15227E+ r 4 

.467862975*01 


39 01 

.377931796+01 

. 12286 7536+Vl 

.399820975- '2 

•142828585+05 

.272553625+ 4 

.4 3 , “ c 69 2*5+ -4 

.17* 7 79555+ <9 

• 1865922 75+ U* 

9 



.122387116+ 71 

.596987145+ C, 

. 1 7 8 C 5 1 6 8 E ♦ 0 2 

. r 

• 7 9 ° 5 74 Q 2 5- 74 

.43957557E+; 4 

.337395156+01 


39 12 

.36956 47 5+.1 

.19 09 7-3265 + - 1 

• 5187232 :E - .'2 

•1316 7845 ♦05 

• 3 c 26 q 54RE+ ‘4 

.44 q 3 ::, 264E+ 0 4 

.164339295+ >9 

.2639091 *6+ }* 

Q 



•12 3375 72 E+ 71 

.597278625+' ? 

.166142125*72 

. 

.1*9855865- -2 

.468227655+ 4 

.3756161:6+01 




A- 2 7 


ITERATION MU«B59 2 DC.-'. OTI'lE = .2*7672 3 2 E - 0 6 TJ! 


V ODE 

X 

Y 

R HO 

U 



G A M 

cv 

T* 1 XY 

39 25 

♦ 3= 1323575+01 

* J 

.246=69595- 2 

. 152729815 +05 



. 12344953E+ 31 

• 5 953 893 4 5+-. 3 

* 

39 37 

.37787 3955+ 01 

•12423 3165*01 

•39r375?15- r 2 

.143166325+05 



.123436375+C 1 

• 5957312 = 5+ . C 

. 1 95 775 065+02 

3949 

.375514535*01 

. 192893445* 9 1 

.5:68 38 225- .2 

.131=' , 7535 + 95 



.123423925*01 

.59 = 92053 5+ .9 

. 1642 A 8245+02 

3961 

.3 8726104 5+01 

. U 

.2* ■'651475- 2 

• 152=22545 + 05 



.123487905+71 

• 59**242 '5 + I 0 

• 

39 73 

.3345354 35+ 11 

.125533555+01 

•382*l?755-:2 

. 143498535+05 



.123474595*01 

.59*761 295+ 00 

.1835 1065= + 02 

39 85 

.383493965+01 

. 19 4959625+ M 

•*959C4=15-C2 

. 1322C2835 + r -5 



. 123467235+ U 

.5949478 8 5+3 0 

. 162185885+82 

3997 

.394224725* Cl 

. ■■> 

•23*6285*5- 32 

. 153657345+C5 



.123526125+01 

. 59345 8 565+ 3 u 

9 

4 0 09 

.391826235+01 

. 126852 915+ 01 

.374847 r- 45-02 

. 143822345 +05 



.123512965+01 

.5=3796-8 35+ 0 6 

. 1 01294285+02 

4 C 2 1 

•390499765*01 

• 197 007245+ Cl 

.*=4094825- 02 

. 13249736^+05 



•123505695+ Cl 

.5 = 39 7*7 05+ CO 

.16 06 16635 **2 

4033 

.4 01215625+01 

. n 

.22 = 2 07685- 32 

.153164865+05 



.123564495+31 

• 3?24922°5 + 0 0 

. 0 

4045 

.398944365*01 

• 128153715+ C 1 

.36 75 49255--' 2 

.14414C205+C5 



.123551475+01 

• c 3 2 8 1 9 6 9 r + . f 

. 99 16 04335+01 

40 5 7 

.39753294 5 +G1 

. 19 9 236 615+ 31 

.*73397375-32 

.132779705+05 



. 123544275+01 

.'93:cC9!5+ r 0 

. 158966025+02 

40 69 

.4 8234785+01 

• - 1 

. 22*627245- 2 

. 15328 T455+ 05 



.123603015+01 

.5 = 1525285+ . 0 

9 i' 

4 C 8 1 

.4 :589 3355+1/1 

. 12 9453875+ 01 

.36' 431825- 2 

. 1 4445 3725 *05 



.12359 9145+01 

.591547 = 55+ 3 C 

.97144 29 85+ Cl 

4093 

.4 ■.4594555 + 01 

.291046055+91 

.*61725645- 32 

. 133 95 7765 + "5 



•12358 3 G3E+ 01 

.592:26395+':5 

. 157593055+02 

41 05 

.415283235+01 

* 0 

• 2 2 9 5 5 u 7 5 - 2 

. 15342.7315+05 



. 123641695+ 01 

.5= 35574 15+ . C 

• L 

4117 

.412965765+01 

.130737675+01 

.353293215-02 

.1*4763=35+95 



. 123628935+ C 1 

. 59 .£ 75155+ O 9 

. 95277 1215+01 

4129 

.4 11 685645+01 

• 2.3.41855+': 1 

. *5122:385- c 2 

. 13332102 r +95 



.123621955+01 

.5 = 1 351 r ’25+ 30 

.156 162485+02 

4141 

.*22362035+01 

. 0 

. 219 3 56985- -3 2 

. 1536 -5435+05 



• 1236S 354 E+ 3 1 

.3 = 95 63555 + . 0 

. c 

4153 

.42 O 0 731 3 E+ 01 

. 13 201 j 3 0 5 + vl 

. 3*61 379 = 5- '2 

. 14 5 071 3 05 +*05 



. 123667 = 8 E+ M 

. 589=014 = 5+ ■. 0 

. =356 37155 + 0-1 

4165 

.4198 072 65 + 01 

• 235C1831E+3-1 

.**1 -19115-02 

.133575 99 5 +05 



♦ 123661035+ 3 1 

.55 '.74695+ C J 

. 15&54 995 5+ 02 

4177 

.429472275+01 

. c 

. 21 5665255-32 

.15379 1855+05 



.123719565+31 

. 5 =6 6 18 5 85+ 3 0 

r 

. 

4169 

.42721 1065+ Cl 

. 13 3271=45+ M 

. 3 3 = 6 91 165- '2 

. 145375885+ 05 



•1237J7155+9 1 

.5 = = 926715+ : 0 

.92 '.78 8145 + Cl 

42 M 

• 42556:535+01 

.2069777:5 + 01 

- *31745 2 35- 3 2 

. 1328 13=45+05 



. 1237 029 5+ 01 

. = =9 0 9725 5+3 -3 

.152? 14495 + "2 

42 1 3 

.4 36615o35+!;l 

. i 

.2134 98 58 5- '2 

• 153951 7 9 r + ? ''5 



. 123758775+ 01 

• r 876473 = 5+ 8 

. 

42 25 

•434391645+ 31 

.134522775+01 

. 351546645-' 2 

. 145677715 + J5 



.123746515+01 

. = = 7=5„715+- 0 

. 9 0766 2 = 45 + "l 


346 7359 .*5- 02 


26 7347725+ 04 


388745975+ 04 


2621 !24or+ 04 


38454987 5 + 04 


256948815+ 04 
0 

38147374 E+ 34 


25 198C835+04 


378147805+ J4 


2&729379E* 04 


37549499 5+ 34 


242971595+04 


372395705+04 


239381195*04 


3693 06855+ 9 4 


2 3 56 6 277 5 + = 4 


3656 '7835+ 04 


SOS 

«U 

•322410735+ '-4 
.74 C5 02 ME- 0 4 
•*3C441535+ 04 

• 7 Q 8965 75E- 0 4 
.4*3860155+ :4 

• 148136 8 65-''- 2 
.3 1 9 8 75= 75+ r -4 
•74C4 qo 275-C4 
.42=847405+ 04 
.798956 *35- 94 
.44=186315+04 
,147 Q 5874E-02 
.315378405+ 04 
.74 r 438 81E- 34 
.429277 HE* 34 
.798 033 06E- C 4 
.4 4 7452 335 + ?4 
.145192395-92 
.312044955*04 

• 74 0486 49E- 04 
.4 2871759E* “4 
.797971°! E- 04 
•446733965+ -4 
.14459396E- 32 
.3 09 .7?? 6 3E* -'4 
.74 0 4 7 8 4 6 E - ?4 
.428154655*04 

• 79 7 048 22E- 04 

• *45944 03 E* 34 
.14 1772=25-02 
.3 06369975+ r 4 
.7* 0475925-94 
.4275745 3 E+ 04 
.7965492 0E- 04 
•*45238255+04 
.140193325-02 
.3 041*8605+ 34 
.740 477515-* 4 
.426965845+ 04 
.7=6324 155- 04 

• 4**6277 0E* '4 
.139348625- 3 2 
.3 021 = 1 76E+ *4 
.740470195- 34 
.4 26321 7 T E+ r 4 
.7948=8915- *4 
.*4 4 0 = 1 c IE* r 4 
. 136M2 62E-0 2 

• 3 '9 05234 = 5 + 94 
.7**492615- 34 
.* 2564 9225+ 34 
.795696765- "4 


.152539795* 09 
•240656325*34 
•170103975*09 
.429 35 32 75* 04 
.164225285+ ^9 
.466911665* ,4 
.151983545*09 
.235531365+04 
. 170146045*09 
.429 '’36275*04 
.164 109085*09 
.4653664 15*0-4 
.151359655*09 
•23232159E* r 4 
.1731 739*5* 09 
.426768645*04 
.164 019685*09 
•463699G4E+G4 
.150739645*09 
.2254 0*46 5*24 
.170221725*09 
.425 5 24175*0 4 
.163911525*39 
.462 "‘6 07 05* 34 
•150204615+ r 9 
•223975715+' 4 
.172263375*39 
.4242746 *5* -'4 
.163818135*59 
.460372465*04 
.1498.08445* ’9 
.217 1359 35 + -4 
.170319 025* 3-9 
.422992685*04 
.16370-6625*09 
.4537=2655*04 
.14954=675*09 
.213 987585* 0 4 
.173 369 155* 29 
.421656255+04 

.163596645+ 39 
,4572 8 791E+'' 4 
.149378275*09 
.2111331*5*64 
.17042031 5+ 3-9 
.420251985*04 
.163472 n 5E+ -9 
.4 55 "58725+ '4 
.14922 9 6 1 5 +' 1 9 
.208 735125+ >4 
• 170470975* .9 
.418776135 + ■- 4 


.646348765+03 
.4737119 55+ Cl 
.182123725+04 
•338354 C6E+C1 
.2571399 .5+04 
.3063686 3E + 01 
.615889265+03 
.479567475+01 
. 1778628 15+04 
. 339356955*01 
.250760335*04 
.307198475+01 
.587193015*03 
•46531 33 4E*ul 
.173824895*04 
,340338675+ Cl 
.243911365*04 
. 3 J8 1 435 65 + 0 1 
.561390135+03 
. 4= 08422 = 5+01 
.169943745*04 
.341311265*01 
. 2375279 3E + G4 
.309 0438 ■■5 + 01 

• 539352"* 75 + 0 3 
.496037795*01 
.166149645*04 

• 34 22 9494 5 + 01 
.230971675*04 

• 3 .'99969 75 + 01 
.521315775*03 
.5308 w56 5 E + C 1 

• 16238 03 75 + 94 
.3432:5765*01 
.2249*20:5*04 
.310864425+01 
.506781735+03 
. 5 C 5 1 0 5 9 E ♦ 0 1 
.158589055*04 
♦34*355795+01 
.219125395+04 
.311693 815 + 01 
. 4947684 25 + 03 
.508921385*01 
.154751655+04 
. 3454518 '-5*: 1 
.213905255+04 
.3123931 25*. - 1 
.*84241665+: 3 
.512283835*31 
. 1 5 r ‘6662 L 5+34 
•346595495* 01 


232728475+ 04 



ITERATION \UM.R5R 20311 


DTIVE 


• 247672'. 2 5-16 


MODE 

X 

Y 

*HG 

J 



GAM 

CY 

T J XY 

42 37 

.433146475+01 

.21*92 0295+ 31 

.423168995- 2 

. 13 4 14 '9 25+15 



.123738735+ .1 

.58 51 186 5+ .!> 

. 15 95 3976 5 + 12 

42 49 

.443791425+01 

• ’1 

•2113971 7 5-C2 

• 15 4 ''5 4615 + 15 



. 123798 155+ 11 

.5866743*5+ ' 9 

• 0 

4261 

•44158 59 85+Ul 

• 135762935+ :i 

.324131275-02 

♦ 1 * 59 7 6995 ♦ "5 



• 12 3 78 6 05 5+ C 1 

. 5 9 6 9 7 3 3 = 5 + • ■• P 

.8963*8585+01 

42 73 

.441366275+01 

.21 9*46335+3! 

• 415393455 - ”2 

. 134254395+05 



•123779355+^1 

.58713*625+ 0 1 

. 148867315+12 

4235 

.451 jR 2535+01 

.0 

.239378455- 02 

. 154 08 7155 + 35 



.123837735+01 

.5857 11825+ 1u 

. 3 

4297 

.44*825225 + 01 

.136992615+01 

•315689725- 02 

.146274125+35 



.123825735+01 

.5359945 8 5 + 1 1 

.386794235+01 

43 09 

.447621045+01 

.212756 075+ -1 

.4383 Jl 1 0,5-0 2 

.134455655+05 



. 123 619175 + 91 

.586157175 + '. 1 

.146656255+02 

4321 

.453249655+ J1 

• 0 

•207631725-12 

.154157825+15 



.123877505+11 

.5847252 95+ 00 

. C 

4333 

.4561 C0525+01 

.138211^45+91 

.309274635-12 

.1*6569345+35 



. 123865715+ 91 

.5859141 85+ “ 9 

.878628495+01 

4345. 

.454911945+01 

.214649765+01 

.4 9,16993 3 5- '2 

.134644255+05 



.123359135+01 

.585174135+ C 0 

. 14435 q l 25+02 

43 57 

.435523835+ J1 

• 0 

.2364231 15- 92 

• I53 Q 917 Q 5+'J5 



. 123917435+ 91 

.582747*65 + ' 

. 1 

*369 

.46341 3045+ Cl 

.139421035+01 

.3119 36 76 5- '2 

. 146862475+05 



.123915845+ 91 

.584-*32175+ 0 9 

.371573365+01 

4381 

.46224 0165+ 01 

•2165276 JE+01 

.395531415- 02 

.134821715+35 



.123899405+11 

. 5 8 4 1 8 9 3 8 5 + 0 9 

.1421^2355+12 

*393 

.472856235+01 

• 9 

.2 959653 15- 0.2 

.15391 9625+75 



.123957665+01 

.582768675+ '9 

. r- 

44 95 

.473762985+31 

»14'62 1' 165+01 

.294725945- 02 

.147152555+05 



. 1239*6 195+ 91 

• 583148125+': 1 

• 865 29 •: 055 ♦ 11 

44 17 

.4696-36875+ 21 

•21828 q 835+ 91 

.399495165- 02 

.1349*9815+05 



.123939835+ ;1 

,5332 0278 5+ 0 0 

. 13971 3’'65 + 02 

4429 

.48 J 216145+31 

. 0 

.2 963 48* 05- -2 

. 153854985+05 



. 123°98 375+ Jl 

.581787515+ 1 1 

« .1 

44 41 

.4731545 75+ 01 

.141818315+01 

.28768 8 355-02 

. 1*7478035+15 



.123936745+11 

.582 1622 0 5+ 0 1 

. 35933 qq 35+ 01 

4453 

.477013325+31 

.22 '236645+ 01 

.383722465-12 

. 1351*7755+05 



. 12393 1435+ 01 

• 5*2214235 + 1 r 

. 1374 '413 r + ->2 

44 65 

.43762 J735+31 

. j 

.2 0 75 08645- 02 

.157799075+05 



. 12403 8715+ 0 1 

.58*3 042*5+ 1 J 

• c 

44 77 

•435536345+01 

.142988675+01 

.28 987** 15- 12 

.1*7716935+05 



. 12*02 7 53 5+ )1 

.58 1 _• 7*1 9 5+ ‘ 9 

.3535 '.0 Q 65 + 01 

44 89 

.43446.735+01 

.222 066245+ H 

•377 q 5~635- . 2 

. I352°*63 r + 05 



.1241 21355+ 1 

.581223575+ \ 1 

.135211405+12 

45 Cl 

.495 3 65405 + Jl 

. 0 

.239282165- 02 

.153733405+15 



.124 0 79555+1 1 

.5798 19755+ . 9 

. 

45 1 2 

• 453959S6E + 01 

. 144158355+.; 1 

.274332* 15- "2 

• 1 * 79 * 7* 75 1 5 



. 12* C6*555+ '1 

. 5 8 " C 8 3 9 2 5 + •. 7 

.34 7271955 +"1 

*525 

.4 = 155 :3 85+ 01 

.223*84825+ ( 1 

.372417755- 2 

• 13542 q 3B5+T5 



. 12*062*65+^1 

• 58 ''2317 1 5 ♦ 9 

. 13?°* 4475 + q 2 

*537 

.5 '.2552295 + 01 

* 

.2114 632 0 5- 92 

. 15363 *-.:5 r + 05 



•12*120655+ H 

.57*. 9 30 3 95+1 9 

. t; 


.346735915-02 


V 

THX2 

36151*935+ n 4 

f> 

r 

22 “'26175E* '4 
356990375+ 04 

r ■ . 

o 

0 

228220745+ 04 
25 1875735+ 14 


2 2 654 2 03 5+ 04 
C 

3*66 17405+ r, 4 

r- 

r, 

2 2 51*521E+ >’4 

ft 

■ 341157255 + "4 

ft 

223937225+ 04 
335849*35+ -04 

2223 1843 5 + 04 
33 '■'46331 E + 14 

22 163727 5 + ?4 
325341675+14 

n 

22 ~* *8455+ 1 4 

3 2 *1 0 272 E* 14 


SOS 

MU 

.4435*6795+14 
.13 2 287 19 E- 0 2 
•2 CQ 031*55+ '■'4 
.7* -’471*65- *4 
.42*923965+ 14 

• 7 Q 5599 495-04 
.443 126 ill 5 + 04 
.133733 46E- 72 
.297731625+ °4 
.74 0487555- )4 
.4 2 4179 '■ C 5 + 14 
.79479749E- 14 
.442769575+ '4 
. 137297 CC E- ^2 

• 296695 1 OE^ 14 
.7 4 1445 2 05- 14 
.423413275+ 34 

• 899966 8 2E- 14 

• 442467 755 + 14 

• 38 6938 7 3E- 32 
.2 Q 6 06 Q .*6E+ 14 
•74142996E-14 
*422635735+ 14 
•895439 56£- 14 
.4*2196625+ 14 
.379352545-12 

• 295975 39 5 + 94 
.741429725-94 
•421855 175* 94 

• 892 944 2°E- 9* 
.44 1937925+ 04 
.^74422625- 32 
.296541 785+ 04 
.74 1435215- 14 
.421 981 565+ 94 
•883799255- 94 
•4416 p 3*95+9* 
•36 c 6 Q 3 a 25-32 
.29779 s3 4 35+ 34 
.741447515- 94 

• *2 n 323 275+ 94 

• 38 56133 a 5- 94 
,441446825+94 
*36 5 2' r 7 , 85“ 32 

• 2 9 ° 7 13 755+ 94 
.741463665- 34 
.4 1 = q 935 15+ ^ 4 
.88 2 295 475- 1 4 
. 441223 T 25+~4 
.36 2 ra 155- ^ 
.3 92 1 48 14 5+ " 1 4 
.74 1431575- 9* 


■ r 
T 

.163341755+09 
.454 7778 25+ f 4 
.149 T. 15*25 + 59 
.2066095 '5 + 1 * 
.173521155+99 
.417234535+ q 4 
.1632021 95+09 

• 453773 0 55 + '9 4 
.1437*2455+39 
.2 04 749855+24 
.17.9572245+09 
.4156395 35+C4 
.163 954795+ '9 
.452 *94665 + 0* 
.148429495+ 09 
.2^3262555+14 
.171625125+09 
.414 026725 + 04 
.162901 325+09 
.452 129995+ n* 
.148 143155+09 
.202339735+04 
.173680135+39 
.412353995+84 
.162738945+09 

• 45 1 4 3 ~ 8 05+94 
.147852765+ 39 
.2 92 1 47345 + 14 
.171736175+39 
.410698545+ .4 
.162575915+08 
.45 '755465+14 
.147998355+ 39 
.212855335+04 
. 179791445+ 99 
.409 059755+14 
.162403535+05 
.45 91 19565+3 4 
.148 ..'13*25+ , ~ q 
.204512565+04 
.170643015+09 
•417453755+04 
. 16222*365+ 09 
.44?45?7?5+ '.a, 
.14*233 3 35+ n 9 
.207768335+94 
,17386 3 775+'' ! 8 
« 4 '5 9 9 59 3 5 + ' 4 

• 162'~'43725+ r 9 
.448865316+04 
.14851 1565+09 
•2l:3* oq 95+G4 


P 

H 

• 2C9111785 + C-4 

. 3 1 3 5 (.• 0 9 7 5 ♦ 0 1 
.47*585925+03 
*515177525+ Cl 
. 14684895E+G4 
.347784325+0 1 
.204814585+14 
.313413365+01 
.465822035+03 
.51753882E+31 
. 143026395+04 
.349 01256E + 01 
.2 009292*5 + 04 
•3133963 q E+Cl 
.458581645+03 
.519246455+01 
•13912833E+94 
.35 0271866 + 11 
. 19734224E+04 
.3142244*5+01 
.4536423 j5 + 0 3 
.520122265+31 
. 1352856 '5 + 04 
.351551495+01 
. 194 016C t'5 + 04 
.314500635+01 

• 45 24 r 4 4 1 5+9 3 
.520026575+01 

• 13152 45 75 + 0 4 
.352638485+01 
. 19 0765265 + 04 
.314761445+91 
•454834875+03 
.518620715+11 
. 1276716 C 5+C4 
. 3541 172 75 + 01 
.187672585+14 

• 314969 8. 65 + 01 
•*6112*975+03 
.5164 5 1935+ Cl 
.124353075+04 
.355372095+91 
.184562425+04 
.315216645+01 
.*77881125+03 
,512951215 + C-l 
•120985135+ 0 4 
.356584155+01 
. 161625015+1* 
.315394425+01 
.483466*55+03 
.508472775+01 




A- 2 9 


ITERATION M U » 5 E 8 2 3 0 C 

" T I u 

E = .2*7672 C2E 

-06 TT M E = .3467359* 

E - 1 2 




NODE 

X 

Y 

PHO 

l’ 

V 

SOS 

r 

P 

I E 



GAM 

cv 

THXY 

THXZ 

MIJ 

T 

M 


45 49 

• 5 'C57671E+31 

• 14 5 3 18 43 5+;l 

. 2 6 = 1 11 25 E - 2 

. 14824 8645 ♦ 05 

• 219. 01 5 335 + )* 

.41839574E+ 24 

.17092914E+ -9 

.117819135+44 




. 1241 J98.3E + M 

.579.91335+ . ■■ 

.84 03 ? 1235 + f 1 

# 

t q 7 ° 59 . 1 9 3 E- ' 4 

• 4 Z 4 4 2 2 5 6 E ♦ ' 4 

. 3577*4675+01 


4561 

•495483565+01 

•2256=6545+ : 1 

.3£6°*?455- 2 

.1 3 c 5 5 364 r + 05 

• 315061 02 5+ '"4 

.4^1 ft 76 74E+ *4 

.161853065+ .9 

.178730825+04 




•124103 8 CE+ Cl 

• 57 92 35 47E + . 0 

. 13 2 846755+ 02 

A 

.356336 12E- 02 

•44926511E+04 

•31 5 5 66365+01 


45 73 

• 51u £85765+01 

• u 

• 2 1 4 0 2 8 5 8 E - ’ 2 

. 15348 0355 + 05 

.0 

.3 05 053275+ 04 

• l *8 8 0127E+ ri 9 

•*°859C12E+03 




•12416 199 5 + *1 

.5772 4 , '555 + 0 9 

• 0 

t r* 

.741502S8E- 04 

• 21439098E+ , " | 4 

• 5 S3 118165 + 01 


4585 

.5 13 1395 15+ni 

• 14646=155+ 3 1 

.262249695-02 

, 1 4 8499205+ 95 

.217319625+04 

.41 8255485+04 

• 1 70964 9 2E + 09 

.114852895+04 

9 



. 12415131E+ U 1 

• 5 7 3 3 9 6 2 4 E ♦ : C 

. 332579 385 + 01 

.0 

.876858*45- 2* 

.403 05243E+04 

. 35 86 26 025+ 01 


4597 

• 5 17 06 15 9 E + Ul 

.227473595+ 01 

• 3616 14*9E- r 2 

. 13 5 666115+05 

,3t 9946415+94 

.4* 0 82770E+ 04 

.16165934E+C9 

.175933255+04 

9 



•124145 3 9E+ Cl 

• 57S23775E + l C 

. 12 = 6 ° 07 = 5 + ' 2 

. 1 

.352 !66 C 7E- 02 

.447751145+04 

,315682615+Cl 


4t J5 

•51766414E+ 01 

. 0 

.216933645- .2 

• 153268 9 25 + C5 

# r 

.3 "8 315 *’2 5+ 04 

.1*9 576535+ 09 

.516032585+03 

1 



. 1242 0358E+ C 1 

.576? 4 75 5 5 + ' C 

• : 


.741532195-14 

.219919695+04 

.4971152 CE + 31 


4621 

*515746415+ Cl 

. 14761 ,47E + 31 

.256772875 - 02 

. 148742485 + 05 

.215328765+04 

.4 1764344F+ 04 

.17 j9 Q 234E+09 

.112C5774E+04 

9 



•124193C5E+ 61 

.577 09553E + 0 0 

. 32365 0 ? * 5 + 1 

.0 

.87*787 1-5- 1* 

.401 73362E+04 

. 359856955+01 


46 33 

.5146855 1E+01 

.229246125+ 01 

.356472455- 02 

.135772835+05 

,3 0 494541 E+ 04 

.** 0623 165+ 04 

.161456445+09 

. 173212255+34 

9 



.12418723E+01 

• 5772 37*15 + 3 0 

. 126 9 85535+ G 2 


• 34 < =167 1 GE- 12 

.447185015+04 

.3158146 15+01 


4645 

•5252893 OE+31 


. 22 3332205- 92 

. 15 3 or 66*5 + B5 

. *» 

.3. 11804 8 15+14 

.149358 a 7E+'-9 

.535943775+03 

1 



.124226825+01 

• 57629417E+ C n 

• >■' 

# 7 

•741563 T 5E- r * 

.223861715* 1 4 

.490712885+01 


*657 

•523461775+01 

•1*8742545+ 01 

.251681415- ,2 

. 14 8 9 = 1 9 ° 5 + ■' 5 

• 212954655+' 4 

.*16954 215+ 0 4 

.171 024695+09 

. 1 09477445+04 




. 124221° 05* vl 

• 57641122F + C 0 

. 31347 47^5 + 91 

. 0 

.872725 14E- *4 

.4 32 076 5+ ‘'4 

.36 C941985 + 01 


4569 

.5223576 25 + 11 

• 231 :3423E+ 31 

.3511 *74*E- -2 

•1358=1195+ 95 

• 3 = 0178515+ 04 

.4* 13 07 195+ 14 

.161 264985+G9 

. 1703362 8 E ♦ U 4 

9 



.124219135+ 31 

. 5 76475° Qr -t- . * 

• 12457277 r +C2 

. 0 

.3*56' 1 730E- 22 

• 446428955* -'4 

.316046135+01 


4681 

.532964135+01 

. 3 

• 224342 0 2E - "2 

• 1527132 CE***? 

r. 

.315633865+14 

.149676 075+39 

.559 097365 + 33 

A 



• 124246 82 E+ 7l 

.5758 1888E+ 0 

. C 

• r 

.7416151 15-14 

.229356725+04 

.46383 02 3E + l1 


4693 

•5311-599E+C1 

• 149865 45 E +r 1 

. 2469783 15- ■ 2 

. 149215255+C5 

.2 1 '’246455+ 04 

.4 16251615+ 1 4 

.171 '.’63355+ 0 9 

• 1 9 70522 15 + 0* 

3 



• 12424 1965+ Cl 

.5759 338=5+ C' r 

.8 02 0 2 6 955 + 01 

r 

.872151365-34 

.398906375+2* 

.36201549E+C1 


4 7 j5 

•53007=355+01 

.232748215+ ;i 

.3*61 0494 E- 02 

. 1359=9465+05 

.29544912 5+ 04 

.43 a 853 785+ "'4 

.161 ■*’74925+0° 

.167517595+34 

T 



.12*239315+01 

.5759975 1E + * 0 

. 122575265+02 

.0 

.345647 -85- 02 

•445437955+04 

.316382195+01 


4717 

.540683575+31 

• J 

.229069135- 02 

. 1524° 8965+ 05 

# * 

.3 19826875+ 04 

•150C6271E+C9 

•586 C 51515+03 

1 



. 12*266955+ 2 1 

.575341275+ 00 

• V 

• c 

.741663955- '4 

.235*52785+04 

.476535315+01 


4729 

.53536 05 15+ 01 

• 15 397923E+ 31 

•242635 9 9E-C 2 

• 14 = 44 0415 + 05 

.2 072 04005+ 14 

.4156 .*6 125+ ■* 4 

.171 10098E+ ,9 

.1 04827C7E+:4 

9 



. 12*262135+01 

.57545423 E + -0 

. 7393= 1755*01 


• 87 ft 990875- ' 4 

.397£>J73?5 + C4 

•363012U8E+01 


4741 

• 537849 5 JE + 31 

.234475 04 5 + 31 

.3405156*5 - .2 

• 136 * = °£3'+ '5 

.2° 14 ’555 5 + 04 

.43 Q 301 **5+ '4 

. 1 6 09943 or* "q 

. 164431745+04 

5 



.12*259555+ 1 

.5755 16725+.- f - 

. 1 2 :83 2375 +02 


.3*3 149655- °2 

.444409375+04 

.316773385+01 


47 53 

•543464555+01 

• IS 

.23 44 63355- . 2 

. 1521 * a 685 + " , 5 

« . 

.324351 115+ 1* 

• 15 0 5 3 = 685 + 29 

.616394675+03 

i 



. 12 428 7 215+ -’1 

.5748612=5+ 05 

r 

.0 

.741728 115- 04 

.242121655+04 

.468966275+01 


4765 

•54S666775+ 11 

•152C8412E* 31 

.238626665- "2 

. 149657765+05 

.2 <*3847215+ 14 

.4 15 01465E+ ?* 

.171 13639E+C9 

. 1 02785235+ 14 




. 1242S2525+ 31 

.574972195+:: 

. 775646.275 + 01 

# r 

.8 7 06°3°25- ' 4 

.396411595+ 4 

.363938175+01 


4777 

.545672525+01 

.236193905+ 01 

.335397565- -'2 

. 1362 '54 =5 + ’5 

.287175335+ *4 

.439880295+04 

• 1 6 1» 7 1 5 7 7 E + 0 9 

.161564455 + '.4 




. 12427?°3E+ 31 

.575 '33545+ ■ * 

.119 ?5 86*5+">2 

« ~ 

• 342926535- 12 

• 44 332326E + C.4 

•31717 057E+vl 


47 89 

•55629 3 575+^1 

• 0 

• 2 * 0 5 1 3 9 7 E - - 2 

. 15 182 . 97=+ '5 

r 

o 

.32°! 244 = 5+ ? 4 

.151 •j°7?2E + C° 

.651413285+03 

i 



. 1243C7615+ 3 1 

.57437=835+ i 

# 1 

+s 

.741791 *05- ' 4 

.2492598 CE + C4 

•4612873 65+01 


4 8 0 1 

•554526295+01 

.15318 : 0 2 5 + C 1 

.23*895965- „■ 2 

• 14 Q 867595 + * , 5 

. 2 " “24281E+04 

•414*63 90E+ 3* 

.171 170985+ ’9 

.100893725+04 




• 1243 33 095 + 01 

.57**97665 + - ° 

.761 0,3 C 525+'01 

• ? 

.86=265365- 0* 

.395 334555+ 04 

•3648C2945+°1 


431? 

.553543895+01 

•237895905+ j! 

.329521945- '.2 

. 136316385+ 05 

.283917975+ 04 

.4 3 8 342 265+ 04 

.16 ".553325+ 59 

. 156319 *'25+04 

•i 



. 12*30 j 4SE+ ul 

.5745 475 7 r + ' L 

.117652885+02 

• J 

.3 38 6 ■'186 5- 0 2 

.442164135+04 

.3176552 0 5 + 01 


48 25 

.564 1771 55+C1 

• 0 

.24698552 5 - '2 

.15154 07 "5 + 05 

r> 

• 

.334 016245+ 04 

.15170=5 >E+09 

. 688 ? 61355 + 03 

X 



•12*328155+ 01 

•573B9332E+ 0 J 

„ 'J 

• 7 

• 7418 52 465- 14 

.2566= 19 * E + ft 4 

.453692615+01 


4337 

• 56244 Z575 + j1 

. 15*267 .>95 + 01 

.2314181 0 E - V 2 

. 15 007 13 05+05 

.l°6433675+04 

.4 1 396" 2 05+ °4 

.171202435+09 

.991415665+33 

9 



• 12 *3 2 3635 + C 1 

. 57*00 (■ 5 = 5+:: 

.74573 .’1=5 + 01 

. ft 

.86 = 9 35*25- ;4 

•39426=3=5+04 

.365615935+0 1 


42 49 

• 56148 51 6 E + 0 1 

.23958*155+ Cl 

• 3 2*357645 - .2 

. 136*16635+ "5 

.2 = CO 63635+ *>4 

.437914 "4E+ •* 4 

. 1 6 ,363d5E+ L-9 

. 15 c 45 c: 855+0 4 

5* 



•12*321125+01 

• 574059655+ . ft 

. 116 01 63 5 5 +0? 

. : 

.337676 * * E - ' 2 

.*41 : - 2 5 8 ::5+ 0* 

. 313 0844 ’ E ♦ 3 1 




A-30 


ITERATION V J ^ B E R 2GG3 

: DTIV 

E = . 247672 Z 2 E 

-56 TIME = .3467359 

I r - 7 2 




MODE 

y 

Y 

RHO 

J 

V 

80S 


P 

15 



G AW 

cv 

T* X Y 

T HXZ 

MJ 

T 

M 


4 3 61 

.572116545+ 31 

• 0 

• 2537 4 5 6 4 5 - ' 2 

. 15 1 26 3 1 4 5 ♦ C 5 

• V* 

.3 38331 865+ 54 

.152332505+09 

.723363925+07 




. 12 4 348 54 5 + Cl 

• 5 734 062 5 E + 0 J 

• U 

• 

.7* 1° 06 91 E- 04 

.254 1 7j5 = E + !'4 

.4*635 = 2 35 + 01 


*3 73 

• 57 0411135+'Jl 

. 155345335+ u 1 

•225133055- : 2 

. 15L-2S716E + 75 

.192515325+0* 

.413*79845+ '4 

.171232745+ *9 

.97*89=595+ 03 

3 



• 124344*0 E+ Jl 

.57351 '"SrE* - 5 

. 73-- 56955E**! 

. r 

.367 39 6915- r 4 

•3=32891?E+‘* 

.3663=1145+01 


43 55 

.569467375+01 

» 24 1 25 8 72 5+ = 1 

•318588 C 35-52 

.l36527 0 =5+"5 

• 2 77 1 5 357 E* 04 

.437262335+ j* 

.1502 J069E+ 09 

.152261785+04 

9 



•124341=45+01 

♦ 5725663 05 + 3 0 

. 1147575*5+12 

* 

.332362705-02 

.439 °4 1355+04 

.313586465+01 


4397 

.58 51 14245+01 

• 0 

.26 059459E-52 

. 15 19S4535 + Q5 

O 

• 3 4 35 85 5 CE + 04 

•152931455+89 

.768303725+03 




• 1243696S5+ 31 

.572915865+00 

• 0 

.0 

.7*1979 C2E- " , * 

.271 5 0 = 295+ =4 

.*39438025+01 


4 = 09 

.578439735+31 

• 1564 1 484 E+ J 1 

• 2 2 5 1 1 C G 6 E - C 2 

. 1 5 04 5 926E + 05 

•188540355+ J 4 

.4 13 C1942E+ 04 

.171 262295+ J9 

.959264855+33 

9 



. 124365325+ 01 

.57"< 5 1846E+ 5 0 

•714249 3 7E+01 

* 0 

.866975 OPE- 0* 

.392347905+04 

.367140025+01 


4921 

• 5 77513545+31 

.24291=725+01 

.3137737 5E- “2 

. 1366 C 5135 + 05 

• 27321671E + 04 

.436779 43E+ 04 

.160 0 02945+' 9 

.1496 0* 615+04 

= 



•12436293E+01 

.57307522E+ 70 

.1131 51925*02 


.33 085° 6 TF- 02 

.438 79 6 * '■5 + f '4 

.316949435+01 


4933 

• 53817 09 75 + 01 

. 0 

•26741368E-02 

.152756555+05 

. * 

.348 03723E+ 04 

.153 46 6355+ 09 

.6 J936D775+03 

1 



. 12439 067E+ 01 

.572422745+ ' 0 

• 5 

. "! 

.742=41 29E- 04 

•278543575+ 04 

.433 '’184P5+91 


49 45 

.5B6527E45+ 31 

• 157475655+ 01 

•2219966 5 E-'2 

•150648455+05 

.1846 14165+14 

.4 1 25709OE+ 54 

.171293435+09 

.9**2 03 6 35+ 03 

9 



. 12438629E+01 

•572523245+ C 0 

.698617835+01 

.* 

,865917 38E- 24 

.3 = 142 = 79E+ 54 

.3676769 LE+01 


4957 

.555619115+01 

•244567215+01 

.3 085 5613 5-v 2 

. 136669215+15 

•26991289E+ C* 

.436304 175+ 04 

.159626275+09 

.146772=15+04 

9 



. 12 438 4 02 E> 01 

•5725753*5+ c- 

• 1117' - 18BE + v2 

. 0 

.327094465- 02 

.437767675+04 

. 3 193333 5 + 0 1 


4369 

•596286715+ 01 

• 0 

.274 0631 3E- ~2 

.15 043 7395 + 05 

• 2 

.352182 105+ 04 

•1539956 15+59 

•649216565+03 

1 



.12441 183E+ Cl 

•5719267*5+ GC 

• v.- 

. ; 

.74 21 "160E- 54 

.285169=95+04 

.427157975+01 


4961 

.5 = 4677l = 5 + =l 

. 1 5 9 5 2 7 f 15+01 

•21 C 06°49E- 52 

.155635925+05 

.16 07 5 029 E+ 04 

.412126715+ 54 

.1713265=5+09 

.92=539825+03 

g 



. 1244=7635+ 01 

•572025145+ C 0 

.68333 0625 + 01 


.864823375-04 

• 3? 052 07 1 5 + 04 

.368612475+01 


4993 

.5937353 4 5 + 01 

•246201265+01 

.3 7 4369755- 02 

. 13674669E+ 05 

•26538980E+ 04 

.435939695+ 0* 

.159612*35+ 09 

.14451*1*5+04 

5 



• 1 2 44 G 5 31 5+ O 1 

• 57 2 9 79 5 75 + 0 0 

.10983 7895+02 

o 

. 

.324225335- 02 

.4369613 05 + 04 

.31=535375+01 


5 3 05 

• 6 -'4 4 691 55 + 31 

. 3 

• 28. 0486335-12 

• 15 01 8 9 75 E + 05 

.4 

• 356=0 41. 0E + 04 

.154471235+09 

.68 7=4 0 J 6 E + 03 

1 



• 124433 14 5 + 01 

.571*27785+ 00 

r< 

• * 

.0 

.742186 01 E- 54 

•2913*2265+04 

.421576475+01 


5 917 

• 6 .2839495 + 01 

• 15957 1345+01 

•216165185- 32 

.151522525+05 

.177164955+ =4 

.4 1 1679 S 15+ 04 

• 1713 6 2575 + ‘=9 

.915120115+03 

9 



.12442=335+=1 

•571524065+1" 

.6687c 61 3E+ 01 

r 

.865178 35E- .* 

.389 607215+ 54 

.36=357375+01 


5 0 29 

• 6 2315375 + 31 

.247S21925+ul 

• 3 0 0 0 3 6 6 9 E - 2 

. 136799585 + 05 

.26 112391 E+ 34 

.435608935+ 04 

.159412395+09 

.1*223=915+0* 

9 



. 12442675E+ Jl 

*5715 7732 E + 0 o 

. 1 08 06 658 E + 02 


.324849 =8 E- 52 

.436223? 0 E+ 54 

.31 = 712125 + :>1 


5141 

.612714=65+ Cl 

• .3 

.2365324 5E- 52 

. 14T97886E+C5 

* 1 

.359483 665+ n 4 

.154923915+0= 

•92*893475+33 

i 



. 124454635+ '1 

.577=25765+ * ? 

. * 

r 

.742214635- '4 

•297'13?35+04 

*4 172 =63 *5 + 0 1 


5C53 

• 6 1 1 1 665 1 E + 01 

. 16 06 3 6 29 E + Ul 

• 2 1 32 78 2 1 E- 02 

. 1512«8S5E+35 

.173559G7E+04 

.4 1 1223 74E+ 0* 

• 1 7 1 * 0 0 2 = 5+ o9 

.90=7428:5+03 

9 



. 12445 359 5 + 01 

.571019=25+ 5 0 

. 65*782245+01 

. 2 

.862462=05- 54 

.338677 585+04 

.3701133 65+31 


50 65 

• 61331 0S4E + G1 

. 249428255+11 

•29662143E-.2 

.136626755+05 

.255817375+04 

.435386555+ 04 

.159183225+09 

•1*0429465+04 

9 



. 124448365+01 

.571572 315+== 

.105° 3 11 95 +72 

# 2 

.318546605- 52 

.4357=2835+04 

.31=71 :4 6 5+31 


5077 

.621025225+01 

. 1 

•29225712E-02 

.149*16725+05 

, ; 

.36259=6 = 5+ 54 

.15537=425+0= 

.95=463075+03 

1 



. 124476255+ C 1 

•57J42C61E+C 9 

9 *- 

• L : 

.742334 1 3E- 04 

• 3 02 1 32 73 E + 04 

*4131739 05 + 01 


5 r - 6 9 

• 61951 004 E * ;!1 

. 161632675+ 01 

• 21 5 3 * 3 4 4 E - 5 2 

. 1513 = 4255+ 05 

.ix 03 1 Q 6 3 E+ 04 

.41 07527=5+ 54 

.17143=175+09 

.956159955+03 

c 



. 124472345+ 0 1 

• 57 =5 12635+ . 

.641882955+51 


.364821735-'-* 

.38771 = 6 ?E+ 04 

. 37 09 02655 + 31 


5101 

.61 3672335+ Jl 

.251 5 2 3 26 5 + 31 

• 29323322E- '2 

• 136847435 + 05 

.2 5 05 7364 E+ 04 

.*352=1 005+ '4 

.15 9=59355 + 3 9 

. 1 386828 15+3* 

9 



.12*47:; 155+11 

.57.0563545 + 30 

. 1 0376 166E + *•2 

. r 

.325277625- 52 

.435255335+04 

.319674315+01 


5113 

.629 4 C4475+ 51 

. J 

.29733' 57E- 52 

. 14 c 7 "9 89 r + =5 

. r- 

• 36530128E+ 54 

.155 8 1315 5+ =9 

•99 J 544635+03 

1 



. 124496 125+01 

•569912215+ '0 



.74232=^65- 5* 

.3 -6-97385+ -4 

.4 0 9 6 2 5 7 = 5 ♦ 0 1 


5125 

.627921935+ 11 

. 16265 5 4 9 E * 51 

.2 5738342E - :;2 

. 151578675+05 

• 1673174 HE* 34 

.*1 0 26 3635+ 5* 

.171*768*5+ '= 

. 2 71 *56465 + 03 

9 



• 1 2 * 4 9 4 2 5 E ♦ 01 

• 5 7 1 u 2 1 0 5 ♦ 0 0 

.6 2=8=91 2" +01 

. r . 

.86 572= 1 55- '4 

.3 = 672 = 565+ -4 

.37171:575+01 


51 37 

.6271 02025+31 

• 25 26 0 3 99 E + il 

.29 52971 IE- '2 

.136=4=^25+05 

.244733375+ 04 

.435081=85+ .4 

.15871727E+0= 

.137194= '5+c* 

a 



. 12*492125+ cl 

.57 , ^ 5 1 825 + ’• 

. 1 0139 2255+02 

. i 

.316117 = 05- ‘2 

, 454 ^ 4 , '53E + - 4 

.3 l = 52 ff l 1 E * u 1 


51 49 

• 6373d 37 j 1 
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•3717 c 638 ^+ -4 

•38 3 721 395+31 


5785 

•792952775+ 31 

•2737431 3 E+ni 

.23 8 52368E- j 2 

• 135228 075 + 95 

•15B85353E^ 34 

.4 339 ■*« 425+ 04 

.153162725+09 

.111731435+ C* 

9 



. 12492439E+ C 1 

.560165325+97 

.66=98773" +01 

• 

.746470765- '2 

.431 1 7 0*55+ ?4 

.3137=4115+01 




A-33 


ITERATION NUMBER 


D T I M 5 


2 A 7 £ 7 2 ;26-06 


T 


5797 

.8043 = 6256 + 01 

59 09 

• 3. 03 53 524 E + Cl 

53 21 

•803152246+ Cl 

53 33 

.814612496+ 01 

53 45 

.813851266+01 

5357 

. 8 1 3 4 3 02 75+ 01 

5369 

• 3249439 35 + C-l 

5381 

.824228375+01 

5393 

.82383024 5+ 01 

5 3 G5 

•835412056+01 

53 1 7 

.634731336+01 

5929 

• 93435564E+01 

5941 

.346302426+01 

5953 

.845363455+ jl 

5365 

• 645 ul 0076+ 31 


GAM 

• V 

.124954126+91 
.18 0329276+ 1 1 

• 12495 2 Cl 6 + 0 1 
. 26 '..059 = 85 + J 1 

• 12495 J 8 5 6 + 0 1 
. 0 

.12498 J 7 2 6 + J 1 
•181167246+ 31 

• 124=78746+ 31 
. 2313613=5+01 
. 124377645+ 01 
. 0 

. 125907656+ .11 
. 18 1995 05 6 + •) 1 

• 12 5 23 5 73 6 + 01 
.282647 036 + Jl 
. 12 5 3 0 4 746 + 31 
. 0 

•125034926+ Cl 
. 16281 2536 + 21 
•125 C 33156+ 01 
•293916616+01 
.125032176+ 11 
® Q 

• 125 06251 6+ C 1 
. 1S36194S6+ 2 1 

• 12 5 06 0856+ r 1 
.28 51698 9 E* 11 


cv 

•21369942E- 02 
•55=4963 ;r + ;.r 
.161227145- -2 

• 5 5 ° 5 4 5 1 0 6 + v 0 
•23614155E-.2 

• 5 5957 , '99E+:C 

• 2C648738E -?2 
.55 9 = 02 C75 + ?? 
.16 0321396-02 

• 55694S45E+ .. ['■ 
.2334136 8 E-'2 

• 5 5 8 9 7 1 i 5 ♦ r ' 
.1995 32525-52 
. 5 5 6 3 9 0 1 5 E + C 0 
. 159*9851 E- 72 
.558342 0 96 + VO 
.231723 0 7 E - v 2 
•558365226+ .0 
. 19 353446E- V2 
•5576=21 ' E + ^ 0 

• 15372751E-';2 
•557731595+ C 

• 23 0 617466 - '2 
.557 7 534 36+ 0 i? 
.1677494 0E-. 0 2 

• 557 776 046+00 
. 1 5794=566- "2 
.557115056+^0 
.229245336-12 


Tm y y 

. 1 5436 62 jE + 05 

» -J 

. 154237726 + 05 
.425 05 6336+ M 
. 135 0° 'j 4 06 + 15 
.549236376+01 
.15*692866+05 
• .. 

.15431 0866+95 
.40652 9226 + 01 
. 1349 = 6736 +95 
. 63271 865E+91 
. 15499 1846 + ' , 5 

. 15433 1 036 + 05 
. 358 '■'79576 + 01 

• 1 34616886 +05 
.611751986+01 
. 15525 •■'3 26 +05 
. 0 

.15445 0616 + 05 
. 369766 116+ 01 
.134323606+05 
. 532141395+01 
. 1 55493396+05 

• 0 

. 1 5452 '0066 + 95 
.352512976+01 
. 134^6 264 r + 05 




*125 059936+ 01 

.557135525+ 00 

.524426355+01 

53 77 

.336724755+Jl 

* 

.1330705:6- 02 

• 1557 13 = 35+05 



.12509 0456+ 01 

.5564577 15+ * 0 

• 0 

5989 

.356127556+01 

.184415615+01 

• 157162715 - '2 

. 154589925 + 05 



. 125:899 06+ 01 

. 556 4 92225 + <-0 

.335846685+21 

6 0 01 

• 955757266 + 01 

.2964 16265+ -.1 

.2276851 :• 6 - v 2 

. 133736265+05 



.12 508 8 04 6 + 01 

.5565113 15+ J 0 

• 53 38 1 1675 + 01 

6013 

.3675929=6+ *1 

*4 

© 

• 17 = 575355 - 0 2 

. 155533925*25 



.125116755+ 01 

.555330=25+:: 

# - 

60 25 

.867 1127375 + 01 

.18 52 0 075 5+ 11 

.1563 06466- 02 

. 15466 7816 + 25 



. 1251 1 73 35+''! 

.55 5? 629 15 +09 

■ 32 095965 e * + 01 

6 0 3 7 

.366721075+ 01 

. 237625625+ v 1 

.226427175- *2 

. 133379345+25 



.125116595+ U 

»55o3SL615+ 

.5267=6=55+01 

6 C 4 = 

• 3 7353 07 7E+ 71 

. 9 

.175357376- 12 

• 156235 23 5 + r 5 



. 12514741^+ -ji 

.5551975 06+ ;C= 

. r - 

6061 

.573 063685+ 01 

.18597*626+ 01 

. 1554 0 C 5 6 6 - 02 

•154733435+05 



. 12514 6 D7E + •- 1 

.55 = 226 = 66 + \ 1 

. 3 268 0 2*75 + 01 

6*73 

• 87778o4 76 + 01 

.288827496+ *1 

• 2 2 51710 ? E - ' 2 

. 1329 = 39 =r + °5 



. 125145 335+ 0 1 

.555243265+* n 

.483?811 a 5+?l 

6 0 85 

.353722515+01 

. T 

• 17233 = 7 *. 5 - '2 

. 15616 1 7 = r 



.125176445+01 

.5545572 75+ 0 0 

• 2 

6097 

.383253315+01 

. 186736=66+ '< 1 

.15*331615- 02 

. 1549 *'69 05 + 05 



.12517 = 225+ 01 

.554584176 + 6 

. 2 = 496 =34 r +01 


6 = .346 7T 5? 76-02 


V 

THXZ 

114633676+04 

153722316+04 
( n 

r- 

1 0966944 6+ 04 
149=75296+94 


l 

1 04726596+ 04 
144280286+34 


998153906+ 03 
1 36948166+ 94 


9518 8 7976 + C 3 
0 

134133016+ 04 


0 

a C 71 3819 6+ 03 
r> 

1242 7 3056+ 94 


.867286726+ ^3 

. 0 

.11828591 6+ 04 


32934595 6+ 03 
11211 791 6+ 04 


79767 = 55 6+ 03 


SOS 

M|J 

• 313 37217 6 + 04 
.74*949836- 94 
.4 02959 376+ 04 
.1 7’ 757596- '3 
.433969756+ 74 
.74 0931 446-02 

• 3 05947Rir + 94 
.744923326- 04 
.4 0298S9C6+ 04 
.1 01936366-03 
*434120406+ j 4 
.737675736- *2 
.298841 18E+ 04 
.744872666- 04 

• 4930413 0E+ 04 
.1=2954856-93 
.434397616+ 94 
.765432316-02 

• 291863 926+ 04 
.75 ''•84= 7 2E- 04 

• 4 03 I 0069E+ 04 
.14 9119246-03 

• 4 3 4 8 29 45 6+ '4 
.169519356-01 
.235 150 = 06+ 04 
•7503==696-'4 
.4 03 15? 345+^4 
. 1 39 018 686 - 03 
.435301 1 1 6 + 04 
.166533336-01 
.278862 076+ 04 
.75 '2 04 6 06- 04 
.4 035 *5 2 06+ 04 

.13 = 397 4Qr.('3 

•43578292E+ 04 
.166421226- 01 
.273114776* 04 
.7= * 917 186- ’4 
.403231356+ n 4 
. 1 3°161 976- 03 
•436320316+ 04 
.166567=86-01 
.268 197 0 35+ 04 
.74=928956- 04 
.403237*96+ .4 
.139384506- *3 
.436 793 296+ 1*4 
.167178 015- 91 
.26* 136 175+ 04 
.74 = 761 835- ■»* 
.403239865+ 04 
.138=90326-03 


T 

•15060526E+ 1 = 
.224 37165E+ ^ 
.171683455+ 09 
.371712 0 8 5 ♦ 0 4 
. 152 0 23 0 75+ '• = 
.431 130595+04 
.1*9 63 0455+ 09 
.214230466+ 34 
.171631696+09 
• 3716889 06+ 04 
.152476446+09 
.431 3375 .-5+ >4 
.14S664295 + r, 9 
.2 04 34 96 1 6+ 04 
. 1 71 68 3 C- 26 + 09 
.371 7051 3E + 04 
.152 r '2CG36+ 59 
.4317=5315+04 

.1948 761 96+ )4 
.171687356+ J= 
.371 7333 15+04 
.151562845+ 0 = 
.432=57365+04 
.146932626+0° 
.1859738 56+ '4 
.171695485+ 59 
.371 75727E+ 7 4 
.151 17255E+T9 
.433 4 02 015+04 
.1*59954 26+ 09 
.1778215 OE + O'4 
•171794356+09 
.371 760325+04 
.150740176+ 09 
.434264366+04 
.1*5232926+ ,= 

. 1 7 !-• 5 2 0 7 2 6 ♦ ■* 
.171714555+09 
.371 724 1 35+ 04 
.15 023 0255+ 09 
.435237 UE + 04 
.1445=0266+09 
.1644 052 0 6 + -4 

• 171 7256 06+ <3 

• 371£5 r 0 ‘ 6 + '4 
•149694445+ 09 
.436 =3 06 55+04 
.144 .'7 0 555+09 
.159*27325+94 
•1717345*5+09 
.371 5125*5+ 04 


P 

M 

• 52 0=93 795 + 03 

.493132936+01 

.651193716+03 

•38381=126+01 

.11 0623415+04 

. 313275925+01 

.480662375+03 

.5656 13475 + 01 

.647*95 3 5 5 + 03 

.383861755+01 

. 1 9 = 395 i; 86 + 04 

.3126*0135+01 

.442933556+03 

.518609365+01 

.644 1998 55 + 03 

. 383920565+01 

.10572*796+04 

.311668055+01 

.403809916+03 

•531927256+C1 

.541134365+03 

.38 39557 96 + 01 

. 108393015+04 

.319513966+01 

» 3 753987 15 + C3 

.5453 .-21 16 + 01 

.636 9331 = 5 + n 3 

• 384 0 90436+01 
.10 7958725 + 0 4 
.3 19424335 + 01 
.353726865+03 
.558354625+01 
•63485996 5+ *3 
.384 9621 =5 + 01 
. 1 07437 J85+G4 
.3 982551 55 + 01 
.33 98906 15+03 

• 57"7634 '5+01 
.631 339656+03 
.384156135+31 
. 1 '7 3826 15 + 04 
.3 06589 = 66 + 01 
.313260025+03 
.5817=25=5+01 
.627555455+ j3 
. 36* 2 785 86 + 01 
.106695156+04 
•305556125+01 
. 298 5475 .5 + 03 
.591216575+01 
.6 232 108 56 + 0 3 
.364445645+01 


13 

1 

9 

= 

1 

5 

0 

1 
9 
9 
1 
9 
a 

1 

9 

9 

1 

9 

3 

1 

9 

9 

1 

9 

9 

1 

= 



A-34 


ITERATION NUMBER 2 0 J ' : D T I w 5 = . 24 76 72 ~ 2 5-96 Tjyr - .346 7359 75- '2 


NODE 

X 

Y 

G A f* 

RHO 

cv 

'J 

T u X Y 

V 

THXZ 

SOS 

u u 

5 

r 

p 

w 

19 

61 09 

.3839945 35+31 

.290.11425+ 01 
. 12517454E+ 31 

.223355465- '2 
.554599 355+ 3 0 

. 13264 226 r ♦C-5 

.458938575+01 

• 1 56474125+94 

.437246425+ 04 
.165913765- 01 

.149206825+09 
.436883945 + -'4 

. 1 06 03:--'2E*It4 
.3 0 4 3 3 3 9 3 £ + 8 1 

3 

6121 

.901 312345+01 

. 0 

. 12520536E+ 91 

• 1 7 n 6 32 3 1 5- 32 

• 553910035 + -.' 9 

.156256575+75 

r. 

• •- 

t 

•2612=1 115+ 04 
.74 = 698 "95- 04 

.143713995+09 

•155974715+^4 

•28 q 1678 a 5+03 

•593C171EE+01 

1 

6133 

.9005876 15+ 01 

• 18 748 746 E+ 01 

* 1252 0 4 75E + 0 1 

•153291?3E-;2 
.55393435E+. 7 

.154877915+05 

.283776845+91 

• 76 771 336 5+73 

.0 

.403165 885+ 04 
.139745945- 03 

.1717373 05 + 69 
.371 3439 05 + 04 

.618526655+03 

.384625965+01 

? 

6145 

.9w035272E+Ul 

•291176995* 01 
.125234145+ 01 

.221384305- '2 
.553947915+0" 

•132312 715+ 05 
.43776 0765 + 01 

.1 C 1288395+ 04 

.437579495+ 04 
.165003745-01 

.148722215+09 

•437446345+04 

•3C325721E+01 

9 

6157 

.912454665+31 

. .) 

•125235675+ 01 

. 169C8836E- 02 
.553255595+1 0 

. 156 3 17.° 25 + 95 

r 

. 0 

.*0 

.259491 775+ 04 
.749655865- 04 

.143482425+09 
• 15379729E+ r -4 

.282571135+03 
• 6023 999 1E + C 1 

1 

6169 

• 912 0 745 05+31 

•138223825+ 01 
• 12523463E+ 0 1 

•1521241 0E- 02 
• 5 5 3 2 7 7 3 1 5 ♦ ’ 0 " 

.154946935+05 

.274723815+01 

.743514815+03 

.403087365+04 
.138 567=3E- 03 

•17173251E+09 
.37111-0.565 + 04 

.613432395+03 
• 38 48426&E + 31 

9 

6181 

.911364255+01 

• 29232271E+ 01- 
.125234135+01 

• 2135165 a E- 3 2 
.55 3-2 893 25 + 3 

. 13201 1255 + 05 
.425125035+01 

.96 13 3273 E+ 03 

.437943 975+ 04 
.163551 25E- 01 

.14«27?83E+w 9 
• 437870.2 ' 5 * v 4 

. 1 0:3966925 + 04 
* 3 0 2 334 8 95 + 01 

9 

6193 

.924 G54Q25+ Cl 

• 'J 

• 12 526 5 9 GE+ C 1 

• 163S4547E-V2 
• 5525 Q 375E + ;'0 

.156336635+05 
• .} 

. * 
• 0 

*25 q 15998E+ 04 
.749653 ME- 04 

.143426915+ '9 
.153367235+ *4 

.28137620 5 + 03 
.603243925+01 

1 

62 05 

*923719035+31 

•13395172E+01 
. 125265 025+31 

• 1 5 0 8 8 5 9 2 E - 0 2 
•55261284E+ 0 0 

. 15501 0045 + 05 
.265955965+01 

.72 9041475+03 

.403009855+ 04 
•13 p 1119 OE-03 

.171 719215+ 99 
• 37 0R779 Q 5 + ‘-’4 

.6 08 058215+03 
• 365 3 45635 + 01 

9 

6217 

•92353377E+ Cl 

• 29 34 51 05 E* 3 1 
. 125264545+31 

.2161 4463E- -2 
• 5526234. 7- 5 + ' 3 

.131697335+05 
.41 081 0145+01 

.94586993 6; + 33 

.43 8 n 41 0 05+04 
.162071525- 01 

.147798625+09 
.433 158195+04 

. 1 029 06 025 + 04 
.301425155+ 31 

G 

6229 

.935315135+01 

. 0 

. 125296545+ ul 

•1688561 ?E-t2 
•55192429E+77 

.156316705+05 

• U 

A 

.0 

• 25°9 a l ”45+ 04 
.749772375-04 

•14350 0915+ C 9 
.154314685+ .>4 

.2831324 05 + 03 
.601238815+01 

1 

6241 

.9355259 55+ 01 

*139664815+91 
• 12529 579E+ u 1 

.14967979E-"2 
• 55 19 4 r 73 E ♦ : 3 

.155 0669 05 + 05 
.25844 C645 + -n 

.699925735+ 03 

. r 

.4 02916665+ 04 
.1 38392465- "3 

.171694395 + 0 9 
.37 0 615455 + "4 

. 6 0 2 7 7 C 6 2 5 ♦ Q 3 
.385252835+01 


6253 

• 9 35366 02 E + 01 

. 29455851 E + ul 
. 12529537E+ - 1 

.212652S5E- 02 
.551 = 49935+ "0 

.131406325*05 
. 405 036405 + 01 

.93 0494125+ 03 

a 

.438194 16E+ 74 
.161416285- "1 

• 147 35555E+ 09 
.433356745+04 

.1 01239475+34 
. 30-0633525+31 


6265 

*947742395+01 

. a 

• 1 2 532 762 E+ 01 

. 17 01 3894E- 02 
•551247C0E+CC 

.1562*6365+05 

,r 

r* 

.262573 145+04 
.749739 125- 04 

.143734615+99 
.157356 015+ 9 4 

.29 79 3585E + C3 
. 595 05844 5 + 61 

1 

62 77 

.9475 Owl 75 + 01 

• 19j3647 0E+‘-'1 
.12532693E+-1 

. 14 64 21 42E - I 2 
.55 126G77E+ 0 

.15511 CJCE+-05 
. 25 O 794405+01 

.679379585+03 

• 4 -02859985+ j 4 
.137150935- 03 

.171 656 8 55+09 
.376418965+04 

.5973861 9i; + C3 
•335391235+01 


6289 

.94736593E+ Jl 

• 29 564 5 45 E + 3 1 
. 12532664E+ jl 

.21 03 38 19 E-: 2 
• 5512693 C 'E + J 0 

. 131 05 3535 + 05 
•38836 °56E +01 

.589684435+ 03 
. .1 

.4 38465 415+ 04 
•157 Q 435 q 5- 01 

.1*68398 9 5 ♦ C 9 
• 43879 n l 9 5+ 04 

•103286 015+34 
. 29°57932E+0 1 

9 

63 n 

• 9 5 9 3 4 2 3 7 E ♦ 3 1 

• V 

•125359 15 E+ 01 

.1716 33515- 72 
•55056166E+ ' " 

.15614 "945 + 05 

r. 

• 

• 0 

.266 183 795+ 09 
.75 C <317995- 04 

.144 188 0 3 5+ -09 
.161672695+04 

.331 81 33 35+0 3 
.586590745+01 

1 

6313 

.959646775+01 

•1910513CE+01 

.12535864E+C1 

.1473948 OE-02 
.55 35 72 7 2 E + 7 

. 15514 04-05 + 05 
.242467365+01 

.656923695+03 

.9 

.4 02831905+ M 
.137894555- 03 

.171605175+09 
.37 0273795 + 014 

.593021635+03 

.385469545+31 

9 

63 25 

•959523595+ 01 

•2967 1 135E+ 01 
.125356365+ j 1 

.2 37 0 Cl 3 IE- 02 
.55 05 78 345+ . <? 

. 13 G71658E+95 
. 380° 47825+01 

.87 039050 5+ 03 

• 43871252E+ 34 
.158325825- 01 

.146359015+09 
.439 1742 = 5+ 04 

.987816885+03 
• 298614565 + 01 


63 37 

.9721 13325+ ul 

• ) 

. 12539114E+G1 

,174547665--:. 2 
.54 8 6 66 03 5 + c 0 

. 155932935 + r -5 

• V 


.271285 155+ "4 
.75 9003265- 04 

.1*4768855+ 0 9 
.1678861 "E+ "4 

. 3 1S4 15465 + 0 3 
.574977775+01 

1 

63 49 

.97 1 9 7104 5 + 01 

. 19172331E+ 1 
• 12539 3 7 6 E ♦ 01 

»14^?36865-:2 
.549976375-*- -0 0 

• 15515438E+05 
.233534875+01 

.632754325+03 

f> 

.402=75425+ n 4 
.136349255-03 

.171545435+ 09 
•37025 c 965+ 04 

.588749615+03 
.3854 3 88 75 *Jl 

3 

6 3 61 

.9713893 35 + 01 

.297755475+ >71 
• 1 2539 054 E* 0 1 

.2053249 IE- 02 
.5498 8 "98 5 + 0 

• 13 0278575 + 05 
.35777 93 3 5 +01 

.81 45535 45 + -0 3 
« r j 

.439234 q 25+ 04 
.154384915-01 

. 145792 09E+ 09 
.44 "1 071 75 +.4 

.951895495+03 
.2971826 -E + C1 


6373 

.9845777 1E+ 01 

. 3 

. 125423615+ >1 

• 177 R 3715E-'12 
.549165 8 “5-*- . ' 

.15579 688 5 + 05 

• > 

m A 

.277 08 7 * 9 r+ :4 
. 7 c -> 4 q 4 215- : 4 

• 145441175 + '' a 

• 175> 99185+0 4 

. 338 354525 + 03 
.562266-15+01 

1 

6385 

.93447*455+01 

. 19238 1 17£ + 0 1 
• 125423 35 EOl 

. 1 4 5 7 2 4 5 4 5 - " 2 
.549 171475-*- '• " 

.15 5 151455+ 05 
.221765535+01 

*6 ' 0 B -3 

• 

.4 0 3 017 62E+ "4 
.137932 345- 03 

. 171*77785+ • ■ 9 
.37* 4241 ,:5+04 

.5-36539575 + 0 3 
.335262915+01 

9 

639 7 

•9 G 4423555+ 01 

.29877715E* 01 
.1254232 5 5+71 

. 2031 34745- 2 
. 5491745 s5 + : "• 

. 12°69 14 = r + 05 
.342687595+01 

.777? 122“E + 03 
• ** 

•439722655+ 04 
.157210215- *1 

•145309155+09 
.4* "97029E+ 04 

.°731 835=5+03 
. 2959232 85+01 

9 

64 '* 9 

.997224715+ 01 

• i 

. 125456565+ J1 

. 181829715- :2 
.549454965+ Vi 

• 15556963 r + 05 

« L 

• 

.28 3688 745* >4 
.7629°1165- 04 

• 14&2 r 9 J0 5 + -9 
.1 83 493375+ J4 

. 362535555+ 03 
•546381415+01 

1 



A- 3 5 


ITERATION NUMBER 2 JO . 

DTI 

ME = .247672.2 

E - 0 6 T 

NODE 

y 

Y 

P HO 

it 



GAM 

cv 

tmxy 

6421 

• 95717«7DE + -n 

. 193C24 HE* 01 

. 145 08559E- :2 

• 15513163E + 35 



. 12545643E+ U 

• 54 3 45776E + .7 

• 2 3=146 55 4 5 + ^1 

64 33 

• 997147 -J 45+01 

. 29977567E+ '>1 

.202 3 810 7 5“ -2 

• 1 29 3 55? IE +85 



•1254563SE+C1 

♦ 5 4 6 4 5 9 2 2 E ♦ :• :j 
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